A micro introduction to Bayesian modeling with BRMS

Julin Maloof

Frequentist vs Bayesian Statistics

Frequentist Statistics

  • What is the probability of the data given a hypothesis
  • Theory relies on data that is not observed
  • Accept or reject null hypothesis
  • Does not require use of priors

Bayesian

  • What is the probability of the hypothesis given the data
  • Does not rely on unobserved data
  • Evaluate ratio of probabilities for different hypotheses
  • Requires specification of prior beliefs (not as bad as it sounds!)

Prior probabilities

  • Bayesian methods require that we specify prior probabilities on the model parameters
  • This might sound like we are cheating
  • Instead we use relatively uninformative priors
  • As long as we have set relatively weak (uninformative) priors, the data will overwhelm the priors
  • (Although if we do have prior information we can use informative priors)

Sampling

  • Bayesian model fitting requires sampling different parameter values
  • For each set of parameters test, calculate the probability of the model, given the data
    • Similar to what we discussed two weeks ago with the nls function
  • Difference is that for nls the goal was just to find the “BEST” set of parameters
  • For Bayesian we need to sample enough to define the probability “landscape”
  • Typically several thousand sets of parameter values are sampled

Posterior probabilities

  • One output of Bayesian model fitting is a set of posterior probabilities
  • These will give us the distribution of possible parameter values for each model parameter, given the data and the priors
  • We can do further tests/comparisons once we have these posterior probabilities

Diagnostics

  • It is critically important to evaluate whether the sampling worked well
  • Typically four sampling “chains” are run
  • Rhat compares variance between and within chains. Should be < 1.05
  • ESS (Effective sampling size) should be > 100 per chain

Let’s try it

library(tidyverse)
library(brms)
dat <- read_csv("../output/height_data_clean.csv")

Please create two objects

  • datsmall that only retains observations from “2023-01-27”
  • datsmallbcd that only retains BH, CC, and DPR observations from “2023-01-27”

Check the number of “cores” on our computer

  • Most computers have multiple “cores” or CPUs that can be used in parallel.
  • Let’s check what we have. We will use this on the next slide.
parallel::detectCores()
[1] 8

What are the default priors?

  • We will use what should now be a familiar formula to specify our model
  • We can ask what priors would be used by default for our formula and data
get_prior(height_cm ~ pop, 
          data = datsmallbcd)
                  prior     class   coef group resp dpar nlpar lb ub
                 (flat)         b                                   
                 (flat)         b  popCC                            
                 (flat)         b popDPR                            
 student_t(3, 3.5, 2.5) Intercept                                   
   student_t(3, 0, 2.5)     sigma                               0   
       source
      default
 (vectorized)
 (vectorized)
      default
      default
  • They “b” class are the coefficients on fixed effects.
  • They have a flat prior (all values from -Infinity to +Infinity are equally likely)
  • We can do better than this while still not biasing the model (next slide)
  • Student_t is like a fat-tailed normal distribution.
    • The intercept is centered at 3.5,
    • sigma (standard deviation) is centered at 0

Student vs normal

Fit a model with height being predicted by pop

m1 <- brm(height_cm ~ pop,
          data = datsmallbcd,
          prior = set_prior("normal(0,10)", class = "b"),
          sample_prior = TRUE,
          cores = 4) # set to 1 if you only have 1 core, or 2 if you have 2.
                       # don't go higher than 4, even if you have more than 4 cores.    
  • We use an uninformative prior for “b”, the coefficient for differences between pops
    • We are setting this prior to be normally distributed, centered on 0, with a standard deviation of 10
    • Note “b” refers to any fixed-effect coefficients. If we had multiple fixed effects we could give them all the same prior, or we could specify them individually.
  • By default brms uses 4 chains cores = 4 means that we want each chain to run in parallel in a different CPU in our computer.
  • By default, each chain has 1,000 warm-up samples and 1,000 real samples

Fit a model with height being predicted by pop

CORRECTED: cores=4 not threads=4

m1 <- brm(height_cm ~ pop,
          data = datsmallbcd,
          prior = set_prior("normal(0,10)", class = "b"),
          sample_prior = TRUE,
          cores = 4)
Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1

The output shows us the progress of each chain

Model summary

summary(m1)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ pop 
   Data: datsmallbcd (Number of observations: 155) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
          Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Intercept     3.22      0.12     2.98     3.46 1.00     3900     2762
popCC         1.46      0.25     0.97     1.95 1.00     4014     2729
popDPR        1.89      0.31     1.30     2.48 1.00     4345     3075

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma     1.25      0.07     1.11     1.40 1.00     4181     2951

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).
  • Estimate is the estimated height
    • Intercept is for the reference population (BH in this case)
    • popCC is the estimated difference between CC and BH
  • 95% confidence intervals show our confidence in the estimate.
    If they do not cross 0, then we are at least 95% confident that our differences are real
  • Be sure to check Rhat and the two ESS stats. What did we want these to be?

Model plots

plot(m1, N=3, ask=FALSE)

  • On the left we see the distribution of the posterior probability
    • Check for unusual pattern, bimodality, etc.
    • If the posterior is near the edge of the range of the prior, consider adjusting the prior
  • On the right we see the value for each parameter for each sample of each chain
    • Should look relatively uniform from left to right and chains should look similar

what are the other priors?

Our posterior plots had two coefficients that we did not specify priors for. What was used?

                  prior     class   coef group resp dpar nlpar lb ub
           normal(0,10)         b                                   
           normal(0,10)         b  popCC                            
           normal(0,10)         b popDPR                            
 student_t(3, 3.5, 2.5) Intercept                                   
   student_t(3, 0, 2.5)     sigma                               0   
       source
         user
 (vectorized)
 (vectorized)
      default
      default

Bayesian hypothesis testing

  • Is popCC really different from popBH?
  • Looking for
    • Post_prob to be near 1 (probability that hypothesis is true)
    • Evid.Ratio > 3 (>3 = “moderate evidence”, > 10 = “strong evidence”)
hypothesis(m1, hypothesis = "popCC > 0")
Hypothesis Tests for class b:
   Hypothesis Estimate Est.Error CI.Lower CI.Upper Evid.Ratio Post.Prob Star
1 (popCC) > 0     1.46      0.25     1.05     1.88        Inf         1    *
---
'CI': 90%-CI for one-sided and 95%-CI for two-sided hypotheses.
'*': For one-sided hypotheses, the posterior probability exceeds 95%;
for two-sided hypotheses, the value tested against lies outside the 95%-CI.
Posterior probabilities of point hypotheses assume equal prior probabilities.

Getting posterior probability for an estimate

Sometime I prefer to do this by “hand”

posterior <- m1 %>% as_draws_df()
dim(posterior)
[1] 4000   12
head(posterior)
# A draws_df: 6 iterations, 1 chains, and 9 variables
  b_Intercept b_popCC b_popDPR sigma prior_Intercept prior_b prior_sigma lprior
1         3.5    0.88      1.8   1.2             2.8   -22.2        2.02   -9.8
2         3.6    1.06      1.4   1.2             1.2    -8.9        0.27   -9.8
3         3.5    1.53      1.7   1.3             6.2   -13.2        5.66   -9.8
4         3.4    0.72      1.7   1.2             6.9    -3.7        5.61   -9.8
5         3.4    1.07      2.0   1.2            -1.7   -15.6        2.55   -9.8
6         3.0    1.83      2.2   1.3             4.0    -1.4        2.05   -9.8
# ... with 1 more variables
# ... hidden reserved variables {'.chain', '.iteration', '.draw'}

How many draws are less than 0?

posterior %>% 
  summarize(CC_greater_than_zero = sum(b_popCC > 0)/n())
# A tibble: 1 × 1
  CC_greater_than_zero
                 <dbl>
1                    1
  • All 4000 posterior samples estimate the the CC coefficient is > 0
  • This is very strong evidence the CC height is greater than BH height

Bayesian hypothesis testing

  • Is popCC the same height as popDPR?
  • Looking for
    • Post_prob to be near 1 (probability that hypothesis is true)
    • Evid.Ratio to > 3 (>3 = “moderate evidence”, > 10 = “strong evidence”)
hypothesis(m1, hypothesis = "popCC = popDPR")
Hypothesis Tests for class b:
            Hypothesis Estimate Est.Error CI.Lower CI.Upper Evid.Ratio
1 (popCC)-(popDPR) = 0    -0.43      0.35    -1.12     0.26      18.44
  Post.Prob Star
1      0.95     
---
'CI': 90%-CI for one-sided and 95%-CI for two-sided hypotheses.
'*': For one-sided hypotheses, the posterior probability exceeds 95%;
for two-sided hypotheses, the value tested against lies outside the 95%-CI.
Posterior probabilities of point hypotheses assume equal prior probabilities.

Fit a model with random effects:

  • Really pop and block should probably both be random, but we will keep pop fixed for now
m2 <- brm(height_cm ~ pop + (1|block),
          data = datsmallbcd,
          prior = set_prior("normal(0,10)", class = "b"),
          sample_prior = TRUE,
          threads = 4)
Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 1).
Chain 1: 
Chain 1: Gradient evaluation took 0.000107 seconds
Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 1.07 seconds.
Chain 1: Adjust your expectations accordingly!
Chain 1: 
Chain 1: 
Chain 1: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 1: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 1: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 1: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 1: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 1: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 1: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 1: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 1: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 1: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 1: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 1: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 1: 
Chain 1:  Elapsed Time: 0.191 seconds (Warm-up)
Chain 1:                0.146 seconds (Sampling)
Chain 1:                0.337 seconds (Total)
Chain 1: 

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 2).
Chain 2: 
Chain 2: Gradient evaluation took 1.1e-05 seconds
Chain 2: 1000 transitions using 10 leapfrog steps per transition would take 0.11 seconds.
Chain 2: Adjust your expectations accordingly!
Chain 2: 
Chain 2: 
Chain 2: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 2: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 2: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 2: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 2: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 2: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 2: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 2: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 2: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 2: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 2: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 2: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 2: 
Chain 2:  Elapsed Time: 0.183 seconds (Warm-up)
Chain 2:                0.153 seconds (Sampling)
Chain 2:                0.336 seconds (Total)
Chain 2: 

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 3).
Chain 3: 
Chain 3: Gradient evaluation took 1.4e-05 seconds
Chain 3: 1000 transitions using 10 leapfrog steps per transition would take 0.14 seconds.
Chain 3: Adjust your expectations accordingly!
Chain 3: 
Chain 3: 
Chain 3: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 3: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 3: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 3: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 3: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 3: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 3: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 3: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 3: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 3: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 3: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 3: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 3: 
Chain 3:  Elapsed Time: 0.2 seconds (Warm-up)
Chain 3:                0.165 seconds (Sampling)
Chain 3:                0.365 seconds (Total)
Chain 3: 

SAMPLING FOR MODEL 'anon_model' NOW (CHAIN 4).
Chain 4: 
Chain 4: Gradient evaluation took 1.2e-05 seconds
Chain 4: 1000 transitions using 10 leapfrog steps per transition would take 0.12 seconds.
Chain 4: Adjust your expectations accordingly!
Chain 4: 
Chain 4: 
Chain 4: Iteration:    1 / 2000 [  0%]  (Warmup)
Chain 4: Iteration:  200 / 2000 [ 10%]  (Warmup)
Chain 4: Iteration:  400 / 2000 [ 20%]  (Warmup)
Chain 4: Iteration:  600 / 2000 [ 30%]  (Warmup)
Chain 4: Iteration:  800 / 2000 [ 40%]  (Warmup)
Chain 4: Iteration: 1000 / 2000 [ 50%]  (Warmup)
Chain 4: Iteration: 1001 / 2000 [ 50%]  (Sampling)
Chain 4: Iteration: 1200 / 2000 [ 60%]  (Sampling)
Chain 4: Iteration: 1400 / 2000 [ 70%]  (Sampling)
Chain 4: Iteration: 1600 / 2000 [ 80%]  (Sampling)
Chain 4: Iteration: 1800 / 2000 [ 90%]  (Sampling)
Chain 4: Iteration: 2000 / 2000 [100%]  (Sampling)
Chain 4: 
Chain 4:  Elapsed Time: 0.191 seconds (Warm-up)
Chain 4:                0.171 seconds (Sampling)
Chain 4:                0.362 seconds (Total)
Chain 4: 

what priors were used for random effects?

prior_summary(m2)
                  prior     class      coef group resp dpar nlpar lb ub
           normal(0,10)         b                                      
           normal(0,10)         b     popCC                            
           normal(0,10)         b    popDPR                            
 student_t(3, 3.5, 2.5) Intercept                                      
   student_t(3, 0, 2.5)        sd                                  0   
   student_t(3, 0, 2.5)        sd           block                  0   
   student_t(3, 0, 2.5)        sd Intercept block                  0   
   student_t(3, 0, 2.5)     sigma                                  0   
       source
         user
 (vectorized)
 (vectorized)
      default
      default
 (vectorized)
 (vectorized)
      default

Check model sampling (do your own)

summary

 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ pop + (1 | block) 
   Data: datsmallbcd (Number of observations: 155) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Group-Level Effects: 
~block (Number of levels: 10) 
              Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(Intercept)     0.50      0.20     0.20     0.95 1.00     1139     1582

Population-Level Effects: 
          Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Intercept     3.23      0.20     2.83     3.63 1.00     1660     1426
popCC         1.42      0.24     0.94     1.89 1.00     4314     2985
popDPR        1.85      0.29     1.29     2.41 1.00     4326     2981

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma     1.18      0.07     1.05     1.33 1.00     3867     3094

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

plots

Is it better than the model without the random effect?

m1 <- add_criterion(m1, "loo")
m2 <- add_criterion(m2, "loo")
loo_compare(m1, m2)
   elpd_diff se_diff
m2  0.0       0.0   
m1 -5.3       3.3   
  • Preferred model is listed first
  • elpd_diff is how much worse a model is relative to the preferred model
  • Want this to be greater (ideally 2X greater) than se_diff to favor the more complex model

Weibull in brms

Data set:

Use all dates for the Weibull modeling

dat <- read_csv("../output/height_data_clean.csv")

dat <- dat %>%
  mutate(day=as.numeric(survey_date-min(survey_date)+30),
         mf=factor(mf))

datbcd <- dat %>% filter(str_detect(pop, "BH|CC|DPR"))

Non-linear model specification in brms

Our formula specification is a bit more complex, so we use the function brmsformula to define the non-linear formula.

I am also going to rename a couple of the parameters to be more clear

f1 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmax + Hmin + k + delta ~ 1,
                  nl=TRUE)
  • The first line should look familiar (Weibull formula)
    • I am dividing k by 100 just to get it on a more reasonable scale
  • The second line is new, and it defines how each parameter should be fit with respect to the predictors in our data set.
  • The notation Hmax + Hmin + k + delta ~ 1 means that we want a single parameter value (Intercept) for each of the four parameters
  • nl=TRUE tells brms that this is a non-linear formula.

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)

Set the priors

What priors should we use?

prior1 <- c(set_prior("normal(5,3)", nlpar="Hmin"),
            set_prior("normal(60,15)",  nlpar="Hmax"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta")
)

Why gamma?

  • Want to constrain k and delta to be > 0.
  • Setting lower bounds causes problems
  • Gamma distribution is positive
tibble(x=seq(0,10,.1),
       gamma_2_2=dgamma(x, 2, 2),
       gamma_20_3=dgamma(x, 20, 3)) %>%
  pivot_longer(-x, values_to = "density") %>%
  ggplot(aes(x=x,y=density, color=name)) +
  geom_line()

Fit the model

fit1 <- brm(formula=f1, data=datbcd, prior=prior1, cores = 4)
Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1

Get the summary

summary(fit1)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmax ~ 1
         Hmin ~ 1
         k ~ 1
         delta ~ 1
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmax_Intercept     65.44      1.48    62.66    68.43 1.00     2040     2171
Hmin_Intercept      4.22      0.54     3.14     5.26 1.00     2690     2731
k_Intercept         0.67      0.01     0.64     0.69 1.01     1800     2008
delta_Intercept     6.09      0.47     5.24     7.05 1.00     1832     2321

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.43      0.26    10.93    11.95 1.00     3576     2775

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

fit1 plots

plot(fit1, N = 3, ask=FALSE)

fit 1 predictions

fit1.prediction <- cbind(datbcd, pred1=predict(fit1)[,"Estimate"])

fit1.prediction %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm), alpha=.1) +
  geom_line(aes(y=pred1), color="blue") 

fit 1 plot predictions per pop

  • Please write code to generate the plot below (you can do this by adding one function call to the code on the previous slide)
  • In your own words, how did it do?

Allow a different Hmax for each pop

  • The plants from different pops grow to different heights. Maybe we should model that.
  • With only 3 pops we will use a fixed effect. (Need 5+ for random effect)
  • See below for how we specify that Hmax should take on a different value for each pop
f2 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + k + delta ~ 1,
                  Hmax ~ pop,
                  nl=TRUE)

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b     popCC                  Hmax       (vectorized)
               (flat)     b    popDPR                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)

Set the priors

  • Do these seem reasonable to you?
prior2 <- c(set_prior("normal(5,3)", nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta")
)

Fit model 2

fit2 <- brm(formula=f2, data=datbcd, prior=prior2, cores = 4)
Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1

Fit 2 summary

Any Problems?

summary(fit2)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         k ~ 1
         delta ~ 1
         Hmax ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      4.01      0.53     2.97     5.02 1.00     2749     2667
k_Intercept         0.66      0.01     0.64     0.68 1.00     2385     2282
delta_Intercept     6.04      0.46     5.20     6.97 1.00     2040     2521
Hmax_Intercept     65.55      1.53    62.70    68.74 1.00     2590     2436
Hmax_popCC         10.99     11.90   -12.75    33.86 1.00     3715     2492
Hmax_popDPR        48.97     11.90    25.40    72.03 1.00     3561     2796

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.27      0.26    10.77    11.79 1.00     4004     2873

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

Fit 2 plots

Any problems?

plot(fit2, N = 3, ask=FALSE)

Fit 2 plot the predictions

datbcd$pred2 <- predict(fit2)[,"Estimate"]

datbcd %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop), alpha=.3) +
  geom_line(aes(y=pred2, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2", guide="none") +
  facet_wrap(~pop)

Fit 2 Extend the CC and DPR curves out for all days

  • It would be nice to get predictions for CC and DPR all the way out to day 200
  • To do this, we create a data frame containing each day and pop for which we would like a prediction
  • The function expand_grid creates a data frame with all combinations of the items it is given (days and pops in this case).
pred.df <- expand_grid(day=min(datbcd$day):max(datbcd$day),
                       pop=unique(datbcd$pop))
pred.df
# A tibble: 537 × 2
     day pop  
   <int> <chr>
 1    30 CC   
 2    30 BH   
 3    30 DPR  
 4    31 CC   
 5    31 BH   
 6    31 DPR  
 7    32 CC   
 8    32 BH   
 9    32 DPR  
10    33 CC   
# ℹ 527 more rows

Fit 2 Extend the CC and DPR curves out for all days (slide 2)

  • Now we use our prediction data frame to generate predictions
fit2.predictions <- cbind(pred.df, prediction=predict(fit2, newdata = pred.df)[,"Estimate"]) %>%
  full_join(datbcd, by=c("day", "pop"))

fit2.predictions
     day pop prediction plantID survey_date block row  col   mf rep height_cm
1     30  CC   3.992822    D14B  2023-01-27    D1   4    B    5  12       3.7
2     30  CC   3.992822    D16B  2023-01-27    D1   6    B    3  10       2.3
3     30  CC   3.992822   D112A  2023-01-27    D1  12    A    7   4       4.1
4     30  CC   3.992822   D228B  2023-01-27    D2  28    B    2  13       4.1
5     30  CC   3.992822   D235A  2023-01-27    D2  35    A    1  11       5.4
6     30  CC   3.992822   D239A  2023-01-27    D2  39    A    5   7       3.6
7     30  CC   3.992822   D239B  2023-01-27    D2  39    B    4   5       5.1
8     30  CC   3.992822    D14C  2023-01-27    D1   4    C    4   9       0.9
9     30  CC   3.992822   F112A  2023-01-27    F1  12    A    6   3       6.5
10    30  CC   3.992822   F120B  2023-01-27    F1  20    B    2   4       6.0
11    30  CC   3.992822   F231A  2023-01-27    F2  31    A    7   8       4.9
12    30  CC   3.992822   F117D  2023-01-27    F1  17    D    3  13       4.7
13    30  CC   3.992822   F224C  2023-01-27    F2  24    C    1   6       3.4
14    30  CC   3.992822   F234D  2023-01-27    F2  34    D    3  12       5.6
15    30  CC   3.992822   H111B  2023-01-27    H1  11    B    7   2       6.8
16    30  CC   3.992822   H115B  2023-01-27    H1  15    B    3   3       5.5
17    30  CC   3.992822   H116A  2023-01-27    H1  16    A    4  10       6.3
18    30  CC   3.992822   H116B  2023-01-27    H1  16    B    6   6       4.4
19    30  CC   3.992822   H221A  2023-01-27    H2  21    A    3   5       4.8
20    30  CC   3.992822   H234B  2023-01-27    H2  34    B    4   7       3.4
21    30  CC   3.992822    H19D  2023-01-27    H1   9    D    2   2       4.8
22    30  CC   3.992822   H236C  2023-01-27    H2  36    C    2   1       5.6
23    30  CC   3.992822    J19D  2023-01-27    J1   9    D    2   6       3.8
24    30  CC   3.992822   J116C  2023-01-27    J1  16    C    6   4       5.6
25    30  CC   3.992822   J236C  2023-01-27    J2  36    C    5   5       4.2
26    30  CC   3.992822   J239C  2023-01-27    J2  39    C    2   3       6.9
27    30  CC   3.992822   J110A  2023-01-27    J1  10    A    3   2       6.3
28    30  CC   3.992822   J112B  2023-01-27    J1  12    B    4   4       4.1
29    30  CC   3.992822   J238B  2023-01-27    J2  38    B    3   9       4.6
30    30  CC   3.992822    L18A  2023-01-27    L1   8    A    4   6       3.4
31    30  CC   3.992822   L232A  2023-01-27    L2  32    A    5   3       5.5
32    30  CC   3.992822   L114D  2023-01-27    L1  14    D    1   1       3.9
33    30  CC   3.992822   L117C  2023-01-27    L1  17    C    7   3       4.2
34    30  BH   4.273811    D15B  2023-01-27    D1   5    B    5  24       3.9
35    30  BH   4.273811    D18B  2023-01-27    D1   8    B    6   7       2.2
36    30  BH   4.273811   D111B  2023-01-27    D1  11    B    7   4       3.7
37    30  BH   4.273811   D119B  2023-01-27    D1  19    B    2   8       3.7
38    30  BH   4.273811   D120A  2023-01-27    D1  20    A    4   2       3.3
39    30  BH   4.273811   D226B  2023-01-27    D2  26    B    5  15       2.7
40    30  BH   4.273811   D230B  2023-01-27    D2  30    B    4   3       3.0
41    30  BH   4.273811   D233B  2023-01-27    D2  33    B    6   8       2.2
42    30  BH   4.273811   D235B  2023-01-27    D2  35    B    5  25       2.6
43    30  BH   4.273811    D16C  2023-01-27    D1   6    C    4  12       3.1
44    30  BH   4.273811    D18D  2023-01-27    D1   8    D    7  17       1.6
45    30  BH   4.273811    D19C  2023-01-27    D1   9    C    1   8       2.1
46    30  BH   4.273811   D117C  2023-01-27    D1  17    C    5  14       3.2
47    30  BH   4.273811   D118D  2023-01-27    D1  18    D    3   5       2.3
48    30  BH   4.273811   D119D  2023-01-27    D1  19    D    5   9       1.4
49    30  BH   4.273811   D222C  2023-01-27    D2  22    C    7  18       1.4
50    30  BH   4.273811   D227C  2023-01-27    D2  27    C    7   5       3.1
51    30  BH   4.273811   D229C  2023-01-27    D2  29    C    4  13       3.3
52    30  BH   4.273811   D229D  2023-01-27    D2  29    D    2   9       2.9
53    30  BH   4.273811   D230D  2023-01-27    D2  30    D    1  10       2.3
54    30  BH   4.273811   D235D  2023-01-27    D2  35    D    3   6       3.7
55    30  BH   4.273811   F110A  2023-01-27    F1  10    A    7  19       2.9
56    30  BH   4.273811   F115A  2023-01-27    F1  15    A    6   9       3.6
57    30  BH   4.273811   F117B  2023-01-27    F1  17    B    2  10       3.4
58    30  BH   4.273811   F227A  2023-01-27    F2  27    A    7   7       3.4
59    30  BH   4.273811   F229B  2023-01-27    F2  29    B    4   5       4.8
60    30  BH   4.273811   F240A  2023-01-27    F2  40    A    3   8       3.6
61    30  BH   4.273811   F240B  2023-01-27    F2  40    B    5  12       6.5
62    30  BH   4.273811    F13C  2023-01-27    F1   3    C    5  16       2.9
63    30  BH   4.273811    F15C  2023-01-27    F1   5    C    5  26       2.9
64    30  BH   4.273811    F17D  2023-01-27    F1   7    D    7   6       3.5
65    30  BH   4.273811   F113C  2023-01-27    F1  13    C    3   7       2.8
66    30  BH   4.273811   F115C  2023-01-27    F1  15    C    4   4       3.4
67    30  BH   4.273811   F120C  2023-01-27    F1  20    C    5  11       5.1
68    30  BH   4.273811   F120D  2023-01-27    F1  20    D    5   1       2.9
69    30  BH   4.273811   F225C  2023-01-27    F2  25    C    7  20       2.4
70    30  BH   4.273811   F231C  2023-01-27    F2  31    C    6  10       3.0
71    30  BH   4.273811   F233C  2023-01-27    F2  33    C    5  17       3.0
72    30  BH   4.273811   F233D  2023-01-27    F2  33    D    2  11       3.6
73    30  BH   4.273811   F235C  2023-01-27    F2  35    C    5  27       3.2
74    30  BH   4.273811   F235D  2023-01-27    F2  35    D    2   1       3.7
75    30  BH   4.273811    H14B  2023-01-27    H1   4    B    2  12       4.5
76    30  BH   4.273811    H18B  2023-01-27    H1   8    B    3   9       3.2
77    30  BH   4.273811   H114A  2023-01-27    H1  14    A    5   3       6.9
78    30  BH   4.273811   H117A  2023-01-27    H1  17    A    4   6       2.3
79    30  BH   4.273811   H118B  2023-01-27    H1  18    B    5  18       3.5
80    30  BH   4.273811   H224B  2023-01-27    H2  24    B    5  29       4.5
81    30  BH   4.273811   H227A  2023-01-27    H2  27    A    5   4       4.5
82    30  BH   4.273811   H229B  2023-01-27    H2  29    B    5  19       3.2
83    30  BH   4.273811   H233A  2023-01-27    H2  33    A    3  10       2.5
84    30  BH   4.273811   H233B  2023-01-27    H2  33    B    4   7       2.9
85    30  BH   4.273811   H239A  2023-01-27    H2  39    A    7   9       3.0
86    30  BH   4.273811   H239B  2023-01-27    H2  39    B    6   2       3.5
87    30  BH   4.273811    H15C  2023-01-27    H1   5    C    5  28       4.6
88    30  BH   4.273811   H110C  2023-01-27    H1  10    C    6  11       4.1
89    30  BH   4.273811   H115D  2023-01-27    H1  15    D    5  13       4.1
90    30  BH   4.273811   H117C  2023-01-27    H1  17    C    2   2       3.5
91    30  BH   4.273811   H119C  2023-01-27    H1  19    C    7  21       4.5
92    30  BH   4.273811   H222D  2023-01-27    H2  22    D    2   3       2.8
93    30  BH   4.273811   H223C  2023-01-27    H2  23    C    7  22       2.8
94    30  BH   4.273811   H227D  2023-01-27    H2  27    D    6  12       2.8
95    30  BH   4.273811   H235C  2023-01-27    H2  35    C    2  13       2.5
96    30  BH   4.273811    J19C  2023-01-27    J1   9    C    5  30       4.4
97    30  BH   4.273811   J110C  2023-01-27    J1  10    C    3   1       1.7
98    30  BH   4.273811   J112C  2023-01-27    J1  12    C    7  23       2.0
99    30  BH   4.273811   J114D  2023-01-27    J1  14    D    4   8       3.3
100   30  BH   4.273811   J116D  2023-01-27    J1  16    D    5   5       2.7
101   30  BH   4.273811   J119C  2023-01-27    J1  19    C    5  20       4.3
102   30  BH   4.273811   J223D  2023-01-27    J2  23    D    4   9       4.3
103   30  BH   4.273811   J225D  2023-01-27    J2  25    D    2   5       1.9
104   30  BH   4.273811   J230C  2023-01-27    J2  30    C    4  14       4.4
105   30  BH   4.273811   J231C  2023-01-27    J2  31    C    5   6       6.1
106   30  BH   4.273811   J231D  2023-01-27    J2  31    D    7  14       2.7
107   30  BH   4.273811   J232D  2023-01-27    J2  32    D    3  12       3.5
108   30  BH   4.273811   J235C  2023-01-27    J2  35    C    7  11       1.5
109   30  BH   4.273811   J236D  2023-01-27    J2  36    D    5  21       3.0
110   30  BH   4.273811    J14B  2023-01-27    J1   4    B    2   4       2.3
111   30  BH   4.273811   J119B  2023-01-27    J1  19    B    3  11       3.8
112   30  BH   4.273811   J120B  2023-01-27    J1  20    B    7  10       2.2
113   30  BH   4.273811   J241B  2023-01-27    J2  41    B    6   4       5.5
114   30  BH   4.273811   J233A  2023-01-27    J2  33    A    3   2       3.0
115   30  BH   4.273811   J237B  2023-01-27    J2  37    B    1   7       3.0
116   30  BH   4.273811    L15B  2023-01-27    L1   5    B    2   6       2.8
117   30  BH   4.273811    L18B  2023-01-27    L1   8    B    7   2       3.3
118   30  BH   4.273811   L110A  2023-01-27    L1  10    A    2  14       1.4
119   30  BH   4.273811   L111B  2023-01-27    L1  11    B    3   3       2.7
120   30  BH   4.273811   L115A  2023-01-27    L1  15    A    5   7       3.0
121   30  BH   4.273811   L120A  2023-01-27    L1  20    A    7  15       2.3
122   30  BH   4.273811    L14D  2023-01-27    L1   4    D    4  10       2.1
123   30  BH   4.273811   L113C  2023-01-27    L1  13    C    3  13       2.4
124   30  BH   4.273811   L113D  2023-01-27    L1  13    D    5  22       3.7
125   30  BH   4.273811   L116C  2023-01-27    L1  16    C    6   5       4.4
126   30  BH   4.273811   L118D  2023-01-27    L1  18    D    7  12       3.7
127   30  BH   4.273811   L228C  2023-01-27    L2  28    C    6   6       3.1
128   30  BH   4.273811   L228D  2023-01-27    L2  28    D    2   7       2.9
129   30  BH   4.273811   L229D  2023-01-27    L2  29    D    3  14       1.2
130   30  BH   4.273811   L230C  2023-01-27    L2  30    C    4  11       4.5
131   30  BH   4.273811   L234C  2023-01-27    L2  34    C    5   8       3.6
132   30  BH   4.273811   L237D  2023-01-27    L2  37    D    7  16       2.6
133   30  BH   4.273811   L240D  2023-01-27    L2  40    D    7  13       3.8
134   30  BH   4.273811   L241D  2023-01-27    L2  41    D    4   1       3.3
135   30  BH   4.273811   L242C  2023-01-27    L2  42    C    7   3       2.0
136   30 DPR   3.827380   D117A  2023-01-27    D1  17    A    6   4       4.1
137   30 DPR   3.827380   D225A  2023-01-27    D2  25    A    2   7       8.5
138   30 DPR   3.827380   D227B  2023-01-27    D2  27    B    4   9       7.0
139   30 DPR   3.827380    F19A  2023-01-27    F1   9    A    2   5       6.5
140   30 DPR   3.827380   F232A  2023-01-27    F2  32    A    6   2       2.5
141   30 DPR   3.827380   F241A  2023-01-27    F2  41    A    2   4       5.5
142   30 DPR   3.827380    F16C  2023-01-27    F1   6    C    6   9       2.5
143   30 DPR   3.827380   F240C  2023-01-27    F2  40    C    7   4       5.4
144   30 DPR   3.827380   H117B  2023-01-27    H1  17    B    5   2       6.4
145   30 DPR   3.827380   H226A  2023-01-27    H2  26    A    4  13       3.5
146   30 DPR   3.827380    H17C  2023-01-27    H1   7    C    2  10       7.8
147   30 DPR   3.827380   H225C  2023-01-27    H2  25    C    2   6       7.5
148   30 DPR   3.827380   J121D  2023-01-27    J1  21    D    4   8       3.9
149   30 DPR   3.827380    J16B  2023-01-27    J1   6    B    6  10       4.1
150   30 DPR   3.827380   J232A  2023-01-27    J2  32    A    2  12       7.5
151   30 DPR   3.827380   J232B  2023-01-27    J2  32    B    6  11       3.2
152   30 DPR   3.827380   L119B  2023-01-27    L1  19    B    4  11       4.3
153   30 DPR   3.827380   L237B  2023-01-27    L2  37    B    6   3       2.3
154   30 DPR   3.827380   L116D  2023-01-27    L1  16    D    6  13       4.5
155   30 DPR   3.827380   L233C  2023-01-27    L2  33    C    2   2       5.4
156   31  CC   4.068987    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
157   31  BH   4.178526    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
158   31 DPR   4.078825    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
159   32  CC   4.142810    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
160   32  BH   3.883001    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
161   32 DPR   3.916884    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
162   33  CC   4.034876    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
163   33  BH   3.724878    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
164   33 DPR   3.703282    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
165   34  CC   3.867184    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
166   34  BH   3.863181    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
167   34 DPR   3.975043    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
168   35  CC   3.883061    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
169   35  BH   4.023801    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
170   35 DPR   4.017210    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
171   36  CC   3.585535    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
172   36  BH   3.980126    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
173   36 DPR   4.063230    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
174   37  CC   4.045533    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
175   37  BH   4.010852    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
176   37 DPR   3.674292    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
177   38  CC   4.049329    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
178   38  BH   4.176141    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
179   38 DPR   4.490353    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
180   39  CC   4.132678    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
181   39  BH   3.960575    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
182   39 DPR   4.326573    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
183   40  CC   3.864127    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
184   40  BH   4.140640    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
185   40 DPR   3.776637    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
186   41  CC   4.319821    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
187   41  BH   3.952420    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
188   41 DPR   4.199894    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
189   42  CC   3.723331    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
190   42  BH   3.720707    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
191   42 DPR   4.266908    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
192   43  CC   3.829944    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
193   43  BH   4.199972    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
194   43 DPR   3.981887    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
195   44  CC   3.866905    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
196   44  BH   4.237729    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
197   44 DPR   4.178732    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
198   45  CC   4.237426    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
199   45  BH   4.029533    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
200   45 DPR   3.877307    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
201   46  CC   4.284264    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
202   46  BH   4.075417    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
203   46 DPR   4.142105    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
204   47  CC   4.361777    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
205   47  BH   3.961679    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
206   47 DPR   4.382056    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
207   48  CC   4.002317    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
208   48  BH   3.757518    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
209   48 DPR   4.443257    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
210   49  CC   3.850574    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
211   49  BH   4.296565    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
212   49 DPR   4.113538    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
213   50  CC   3.991818    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
214   50  BH   3.986769    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
215   50 DPR   3.948452    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
216   51  CC   4.284538    D14B  2023-02-17    D1   4    B    5  12       3.4
217   51  CC   4.284538   D112A  2023-02-17    D1  12    A    7   4       4.3
218   51  CC   4.284538   D228B  2023-02-17    D2  28    B    2  13       4.0
219   51  CC   4.284538   D235A  2023-02-17    D2  35    A    1  11       5.5
220   51  CC   4.284538   D239A  2023-02-17    D2  39    A    5   7       3.8
221   51  CC   4.284538   D239B  2023-02-17    D2  39    B    4   5       5.2
222   51  CC   4.284538   F112A  2023-02-17    F1  12    A    6   3       6.1
223   51  CC   4.284538   F120B  2023-02-17    F1  20    B    2   4       5.6
224   51  CC   4.284538   F231A  2023-02-17    F2  31    A    7   8       4.1
225   51  CC   4.284538   F117D  2023-02-17    F1  17    D    3  13       4.6
226   51  CC   4.284538   F234D  2023-02-17    F2  34    D    3  12       5.5
227   51  CC   4.284538   H111B  2023-02-17    H1  11    B    7   2       7.6
228   51  CC   4.284538   H115B  2023-02-17    H1  15    B    3   3       6.1
229   51  CC   4.284538   H116A  2023-02-17    H1  16    A    4  10       6.1
230   51  CC   4.284538   H234B  2023-02-17    H2  34    B    4   7       2.9
231   51  CC   4.284538    H19D  2023-02-17    H1   9    D    2   2       4.0
232   51  CC   4.284538   H236C  2023-02-17    H2  36    C    2   1       6.2
233   51  CC   4.284538   J110A  2023-02-17    J1  10    A    3   2       6.5
234   51  CC   4.284538   J112B  2023-02-17    J1  12    B    4   4       4.5
235   51  CC   4.284538   J228A  2023-02-17    J2  28    A    5   6       3.6
236   51  CC   4.284538   J238B  2023-02-17    J2  38    B    3   9       3.5
237   51  CC   4.284538    J19D  2023-02-17    J1   9    D    2   6       3.1
238   51  CC   4.284538   J236C  2023-02-17    J2  36    C    5   5       4.9
239   51  CC   4.284538   J239C  2023-02-17    J2  39    C    2   3       7.7
240   51  CC   4.284538    L18A  2023-02-17    L1   8    A    4   6       4.1
241   51  CC   4.284538   L232A  2023-02-17    L2  32    A    5   3       5.9
242   51  CC   4.284538   L117C  2023-02-17    L1  17    C    7   3       4.9
243   51  BH   4.315813    D15B  2023-02-17    D1   5    B    5  24       4.3
244   51  BH   4.315813    D18B  2023-02-17    D1   8    B    6   7       2.9
245   51  BH   4.315813   D111B  2023-02-17    D1  11    B    7   4       3.5
246   51  BH   4.315813   D119B  2023-02-17    D1  19    B    2   8       3.3
247   51  BH   4.315813   D120A  2023-02-17    D1  20    A    4   2       3.9
248   51  BH   4.315813   D226B  2023-02-17    D2  26    B    5  15       3.0
249   51  BH   4.315813   D230B  2023-02-17    D2  30    B    4   3       3.9
250   51  BH   4.315813   D233B  2023-02-17    D2  33    B    6   8       2.2
251   51  BH   4.315813   D235B  2023-02-17    D2  35    B    5  25       1.9
252   51  BH   4.315813    D16C  2023-02-17    D1   6    C    4  12       4.7
253   51  BH   4.315813    D18D  2023-02-17    D1   8    D    7  17       1.5
254   51  BH   4.315813    D19C  2023-02-17    D1   9    C    1   8      24.0
255   51  BH   4.315813   D117C  2023-02-17    D1  17    C    5  14       2.9
256   51  BH   4.315813   D118D  2023-02-17    D1  18    D    3   5       2.2
257   51  BH   4.315813   D119D  2023-02-17    D1  19    D    5   9       0.9
258   51  BH   4.315813   D227C  2023-02-17    D2  27    C    7   5       3.4
259   51  BH   4.315813   D229C  2023-02-17    D2  29    C    4  13       3.4
260   51  BH   4.315813   D229D  2023-02-17    D2  29    D    2   9       3.3
261   51  BH   4.315813   D230D  2023-02-17    D2  30    D    1  10       2.1
262   51  BH   4.315813   D235D  2023-02-17    D2  35    D    3   6       4.0
263   51  BH   4.315813   F115A  2023-02-17    F1  15    A    6   9       4.1
264   51  BH   4.315813   F117B  2023-02-17    F1  17    B    2  10       3.0
265   51  BH   4.315813   F229B  2023-02-17    F2  29    B    4   5       5.1
266   51  BH   4.315813   F240B  2023-02-17    F2  40    B    5  12       4.2
267   51  BH   4.315813    F15C  2023-02-17    F1   5    C    5  26       2.9
268   51  BH   4.315813    F17D  2023-02-17    F1   7    D    7   6       3.1
269   51  BH   4.315813   F113C  2023-02-17    F1  13    C    3   7       2.9
270   51  BH   4.315813   F115C  2023-02-17    F1  15    C    4   4       3.5
271   51  BH   4.315813   F120D  2023-02-17    F1  20    D    5   1       1.7
272   51  BH   4.315813   F225C  2023-02-17    F2  25    C    7  20       2.2
273   51  BH   4.315813   F233C  2023-02-17    F2  33    C    5  17       3.1
274   51  BH   4.315813   F233D  2023-02-17    F2  33    D    2  11       3.7
275   51  BH   4.315813   F235C  2023-02-17    F2  35    C    5  27       3.5
276   51  BH   4.315813   F235D  2023-02-17    F2  35    D    2   1       3.6
277   51  BH   4.315813    H14B  2023-02-17    H1   4    B    2  12       5.2
278   51  BH   4.315813    H18B  2023-02-17    H1   8    B    3   9       3.3
279   51  BH   4.315813   H114A  2023-02-17    H1  14    A    5   3       7.4
280   51  BH   4.315813   H117A  2023-02-17    H1  17    A    4   6       3.5
281   51  BH   4.315813   H118B  2023-02-17    H1  18    B    5  18       3.2
282   51  BH   4.315813   H224B  2023-02-17    H2  24    B    5  29       4.6
283   51  BH   4.315813   H227A  2023-02-17    H2  27    A    5   4       3.9
284   51  BH   4.315813   H229B  2023-02-17    H2  29    B    5  19       3.1
285   51  BH   4.315813   H233B  2023-02-17    H2  33    B    4   7       2.5
286   51  BH   4.315813   H239B  2023-02-17    H2  39    B    6   2       3.9
287   51  BH   4.315813    H15C  2023-02-17    H1   5    C    5  28       4.7
288   51  BH   4.315813   H110C  2023-02-17    H1  10    C    6  11       4.0
289   51  BH   4.315813   H115D  2023-02-17    H1  15    D    5  13       4.5
290   51  BH   4.315813   H119C  2023-02-17    H1  19    C    7  21       4.1
291   51  BH   4.315813   H222D  2023-02-17    H2  22    D    2   3       2.4
292   51  BH   4.315813   H223C  2023-02-17    H2  23    C    7  22       3.2
293   51  BH   4.315813   H227D  2023-02-17    H2  27    D    6  12       3.8
294   51  BH   4.315813   H235C  2023-02-17    H2  35    C    2  13       2.2
295   51  BH   4.315813    J14B  2023-02-17    J1   4    B    2   4       2.6
296   51  BH   4.315813   J119B  2023-02-17    J1  19    B    3  11       4.2
297   51  BH   4.315813   J233A  2023-02-17    J2  33    A    3   2       3.1
298   51  BH   4.315813   J237B  2023-02-17    J2  37    B    1   7       2.7
299   51  BH   4.315813   J241B  2023-02-17    J2  41    B    6   4       4.5
300   51  BH   4.315813    J19C  2023-02-17    J1   9    C    5  30       5.5
301   51  BH   4.315813   J110C  2023-02-17    J1  10    C    3   1       1.9
302   51  BH   4.315813   J114D  2023-02-17    J1  14    D    4   8       3.4
303   51  BH   4.315813   J119C  2023-02-17    J1  19    C    5  20       5.3
304   51  BH   4.315813   J223D  2023-02-17    J2  23    D    4   9       1.7
305   51  BH   4.315813   J225D  2023-02-17    J2  25    D    2   5       1.9
306   51  BH   4.315813   J230C  2023-02-17    J2  30    C    4  14       5.3
307   51  BH   4.315813   J231C  2023-02-17    J2  31    C    5   6       7.1
308   51  BH   4.315813   J231D  2023-02-17    J2  31    D    7  14       4.8
309   51  BH   4.315813   J232D  2023-02-17    J2  32    D    3  12       3.4
310   51  BH   4.315813   J236D  2023-02-17    J2  36    D    5  21      25.0
311   51  BH   4.315813    L15B  2023-02-17    L1   5    B    2   6       2.9
312   51  BH   4.315813    L18B  2023-02-17    L1   8    B    7   2       2.7
313   51  BH   4.315813   L110A  2023-02-17    L1  10    A    2  14       0.6
314   51  BH   4.315813   L111B  2023-02-17    L1  11    B    3   3       2.1
315   51  BH   4.315813   L120A  2023-02-17    L1  20    A    7  15       2.3
316   51  BH   4.315813    L14D  2023-02-17    L1   4    D    4  10       2.2
317   51  BH   4.315813   L113C  2023-02-17    L1  13    C    3  13       2.2
318   51  BH   4.315813   L113D  2023-02-17    L1  13    D    5  22       4.7
319   51  BH   4.315813   L116C  2023-02-17    L1  16    C    6   5       4.4
320   51  BH   4.315813   L228C  2023-02-17    L2  28    C    6   6       2.4
321   51  BH   4.315813   L228D  2023-02-17    L2  28    D    2   7       2.8
322   51  BH   4.315813   L229D  2023-02-17    L2  29    D    3  14       1.9
323   51  BH   4.315813   L230C  2023-02-17    L2  30    C    4  11       4.4
324   51  BH   4.315813   L234C  2023-02-17    L2  34    C    5   8       4.1
325   51  BH   4.315813   L237D  2023-02-17    L2  37    D    7  16       3.6
326   51  BH   4.315813   L240D  2023-02-17    L2  40    D    7  13       4.1
327   51  BH   4.315813   L241D  2023-02-17    L2  41    D    4   1       3.5
328   51 DPR   4.183179   D225A  2023-02-17    D2  25    A    2   7       8.9
329   51 DPR   4.183179   D227B  2023-02-17    D2  27    B    4   9       7.2
330   51 DPR   4.183179    F19A  2023-02-17    F1   9    A    2   5       6.0
331   51 DPR   4.183179   F232A  2023-02-17    F2  32    A    6   2       2.0
332   51 DPR   4.183179    F16C  2023-02-17    F1   6    C    6   9       2.6
333   51 DPR   4.183179   F240C  2023-02-17    F2  40    C    7   4       6.3
334   51 DPR   4.183179   H117B  2023-02-17    H1  17    B    5   2       8.1
335   51 DPR   4.183179    H17C  2023-02-17    H1   7    C    2  10       8.2
336   51 DPR   4.183179   H225C  2023-02-17    H2  25    C    2   6       7.8
337   51 DPR   4.183179    J16B  2023-02-17    J1   6    B    6  10       3.6
338   51 DPR   4.183179   J232B  2023-02-17    J2  32    B    6  11       2.1
339   51 DPR   4.183179   L119B  2023-02-17    L1  19    B    4  11       3.7
340   51 DPR   4.183179   L237B  2023-02-17    L2  37    B    6   3       4.0
341   51 DPR   4.183179   L116D  2023-02-17    L1  16    D    6  13       4.0
342   51 DPR   4.183179   L233C  2023-02-17    L2  33    C    2   2       5.7
343   52  CC   4.053981    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
344   52  BH   4.217283    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
345   52 DPR   4.012661    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
346   53  CC   4.185297    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
347   53  BH   4.503634    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
348   53 DPR   4.321795    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
349   54  CC   4.159477    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
350   54  BH   4.254655    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
351   54 DPR   4.145688    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
352   55  CC   4.222352    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
353   55  BH   4.194888    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
354   55 DPR   4.409792    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
355   56  CC   4.150925    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
356   56  BH   4.066326    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
357   56 DPR   4.167857    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
358   57  CC   4.488152    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
359   57  BH   4.255801    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
360   57 DPR   4.360911    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
361   58  CC   4.235308    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
362   58  BH   4.246021    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
363   58 DPR   4.166419    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
364   59  CC   4.332433    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
365   59  BH   4.197437    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
366   59 DPR   4.581395    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
367   60  CC   4.392791    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
368   60  BH   4.426652    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
369   60 DPR   4.376965    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
370   61  CC   4.283974    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
371   61  BH   4.651250    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
372   61 DPR   4.328087    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
373   62  CC   4.749299    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
374   62  BH   4.166820    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
375   62 DPR   4.835346    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
376   63  CC   4.261104    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
377   63  BH   4.137865    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
378   63 DPR   4.538657    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
379   64  CC   4.154024    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
380   64  BH   4.158154    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
381   64 DPR   4.607066    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
382   65  CC   4.601989   D112A  2023-03-03    D1  12    A    7   4       5.7
383   65  CC   4.601989   D235A  2023-03-03    D2  35    A    1  11       5.9
384   65  CC   4.601989   D239A  2023-03-03    D2  39    A    5   7       3.9
385   65  CC   4.601989   D239B  2023-03-03    D2  39    B    4   5       5.8
386   65  CC   4.601989   F231A  2023-03-03    F2  31    A    7   8       5.4
387   65  CC   4.601989   F117D  2023-03-03    F1  17    D    3  13       5.7
388   65  CC   4.601989   H111B  2023-03-03    H1  11    B    7   2       8.6
389   65  CC   4.601989   H115B  2023-03-03    H1  15    B    3   3       6.4
390   65  CC   4.601989   H116A  2023-03-03    H1  16    A    4  10       8.0
391   65  CC   4.601989   H234B  2023-03-03    H2  34    B    4   7       3.5
392   65  CC   4.601989   H236C  2023-03-03    H2  36    C    2   1       5.9
393   65  CC   4.601989   J112B  2023-03-03    J1  12    B    4   4       5.8
394   65  CC   4.601989   J228A  2023-03-03    J2  28    A    5   6       4.3
395   65  CC   4.601989   J236C  2023-03-03    J2  36    C    5   5       5.2
396   65  CC   4.601989   J239C  2023-03-03    J2  39    C    2   3       8.2
397   65  CC   4.601989    L18A  2023-03-03    L1   8    A    4   6       5.5
398   65  CC   4.601989   L232A  2023-03-03    L2  32    A    5   3       6.3
399   65  CC   4.601989   L117C  2023-03-03    L1  17    C    7   3       4.8
400   65  BH   4.394426    D18B  2023-03-03    D1   8    B    6   7       3.3
401   65  BH   4.394426   D120A  2023-03-03    D1  20    A    4   2       6.5
402   65  BH   4.394426   D226B  2023-03-03    D2  26    B    5  15       4.2
403   65  BH   4.394426   D230B  2023-03-03    D2  30    B    4   3       4.9
404   65  BH   4.394426   D233B  2023-03-03    D2  33    B    6   8       4.0
405   65  BH   4.394426    D16C  2023-03-03    D1   6    C    4  12       5.9
406   65  BH   4.394426    D19C  2023-03-03    D1   9    C    1   8       2.7
407   65  BH   4.394426   D119D  2023-03-03    D1  19    D    5   9       1.4
408   65  BH   4.394426   D229C  2023-03-03    D2  29    C    4  13       4.7
409   65  BH   4.394426   D229D  2023-03-03    D2  29    D    2   9       4.8
410   65  BH   4.394426   D235D  2023-03-03    D2  35    D    3   6       5.2
411   65  BH   4.394426   F115A  2023-03-03    F1  15    A    6   9       5.5
412   65  BH   4.394426   F229B  2023-03-03    F2  29    B    4   5       5.0
413   65  BH   4.394426   F240B  2023-03-03    F2  40    B    5  12       4.3
414   65  BH   4.394426    F15C  2023-03-03    F1   5    C    5  26       2.9
415   65  BH   4.394426   F113C  2023-03-03    F1  13    C    3   7       3.8
416   65  BH   4.394426   F115C  2023-03-03    F1  15    C    4   4       4.9
417   65  BH   4.394426   F235C  2023-03-03    F2  35    C    5  27       4.2
418   65  BH   4.394426   F235D  2023-03-03    F2  35    D    2   1       4.5
419   65  BH   4.394426    H14B  2023-03-03    H1   4    B    2  12       5.1
420   65  BH   4.394426    H18B  2023-03-03    H1   8    B    3   9       3.5
421   65  BH   4.394426   H114A  2023-03-03    H1  14    A    5   3       6.9
422   65  BH   4.394426   H117A  2023-03-03    H1  17    A    4   6       3.3
423   65  BH   4.394426   H118B  2023-03-03    H1  18    B    5  18       5.0
424   65  BH   4.394426   H224B  2023-03-03    H2  24    B    5  29       9.3
425   65  BH   4.394426   H227A  2023-03-03    H2  27    A    5   4       3.9
426   65  BH   4.394426   H229B  2023-03-03    H2  29    B    5  19       4.4
427   65  BH   4.394426   H233B  2023-03-03    H2  33    B    4   7       2.5
428   65  BH   4.394426   H239B  2023-03-03    H2  39    B    6   2       3.5
429   65  BH   4.394426    H15C  2023-03-03    H1   5    C    5  28       5.3
430   65  BH   4.394426   H110C  2023-03-03    H1  10    C    6  11       3.7
431   65  BH   4.394426   H115D  2023-03-03    H1  15    D    5  13       4.8
432   65  BH   4.394426   H119C  2023-03-03    H1  19    C    7  21       4.6
433   65  BH   4.394426   H222D  2023-03-03    H2  22    D    2   3       3.1
434   65  BH   4.394426   H223C  2023-03-03    H2  23    C    7  22       3.4
435   65  BH   4.394426   H235C  2023-03-03    H2  35    C    2  13       2.8
436   65  BH   4.394426    J14B  2023-03-03    J1   4    B    2   4       3.6
437   65  BH   4.394426   J119B  2023-03-03    J1  19    B    3  11       5.6
438   65  BH   4.394426   J233A  2023-03-03    J2  33    A    3   2       4.6
439   65  BH   4.394426   J237B  2023-03-03    J2  37    B    1   7       4.8
440   65  BH   4.394426    J19C  2023-03-03    J1   9    C    5  30       5.6
441   65  BH   4.394426   J110C  2023-03-03    J1  10    C    3   1       2.2
442   65  BH   4.394426   J114D  2023-03-03    J1  14    D    4   8       3.5
443   65  BH   4.394426   J119C  2023-03-03    J1  19    C    5  20       5.2
444   65  BH   4.394426   J223D  2023-03-03    J2  23    D    4   9       4.9
445   65  BH   4.394426   J225D  2023-03-03    J2  25    D    2   5       1.9
446   65  BH   4.394426   J230C  2023-03-03    J2  30    C    4  14       6.3
447   65  BH   4.394426   J231C  2023-03-03    J2  31    C    5   6       7.8
448   65  BH   4.394426   J231D  2023-03-03    J2  31    D    7  14       7.7
449   65  BH   4.394426   J232D  2023-03-03    J2  32    D    3  12       3.5
450   65  BH   4.394426    L15B  2023-03-03    L1   5    B    2   6       3.8
451   65  BH   4.394426    L18B  2023-03-03    L1   8    B    7   2       2.3
452   65  BH   4.394426   L110A  2023-03-03    L1  10    A    2  14       1.2
453   65  BH   4.394426   L111B  2023-03-03    L1  11    B    3   3       2.5
454   65  BH   4.394426    L14D  2023-03-03    L1   4    D    4  10       3.6
455   65  BH   4.394426   L113C  2023-03-03    L1  13    C    3  13       2.4
456   65  BH   4.394426   L113D  2023-03-03    L1  13    D    5  22       4.6
457   65  BH   4.394426   L116C  2023-03-03    L1  16    C    6   5       4.8
458   65  BH   4.394426   L228C  2023-03-03    L2  28    C    6   6       2.2
459   65  BH   4.394426   L228D  2023-03-03    L2  28    D    2   7       3.5
460   65  BH   4.394426   L229D  2023-03-03    L2  29    D    3  14       2.1
461   65  BH   4.394426   L230C  2023-03-03    L2  30    C    4  11       5.1
462   65  BH   4.394426   L234C  2023-03-03    L2  34    C    5   8       4.4
463   65  BH   4.394426   L240D  2023-03-03    L2  40    D    7  13       4.3
464   65  BH   4.394426   L241D  2023-03-03    L2  41    D    4   1       3.2
465   65 DPR   4.736533   D225A  2023-03-03    D2  25    A    2   7       9.3
466   65 DPR   4.736533   D227B  2023-03-03    D2  27    B    4   9       7.9
467   65 DPR   4.736533   F232A  2023-03-03    F2  32    A    6   2       2.7
468   65 DPR   4.736533    F16C  2023-03-03    F1   6    C    6   9       4.3
469   65 DPR   4.736533   F240C  2023-03-03    F2  40    C    7   4       5.8
470   65 DPR   4.736533   H117B  2023-03-03    H1  17    B    5   2      10.1
471   65 DPR   4.736533    H17C  2023-03-03    H1   7    C    2  10       9.2
472   65 DPR   4.736533   H225C  2023-03-03    H2  25    C    2   6       7.6
473   65 DPR   4.736533    J16B  2023-03-03    J1   6    B    6  10       4.5
474   65 DPR   4.736533   J232B  2023-03-03    J2  32    B    6  11       4.4
475   65 DPR   4.736533   L237B  2023-03-03    L2  37    B    6   3       5.1
476   66  CC   4.660378    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
477   66  BH   4.384649    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
478   66 DPR   4.879122    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
479   67  CC   4.277645    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
480   67  BH   4.488445    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
481   67 DPR   4.811689    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
482   68  CC   4.651261    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
483   68  BH   4.447075    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
484   68 DPR   4.790984    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
485   69  CC   4.739512    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
486   69  BH   4.300966    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
487   69 DPR   4.844630    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
488   70  CC   4.554620    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
489   70  BH   4.711490    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
490   70 DPR   4.854102    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
491   71  CC   4.995747    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
492   71  BH   4.681415    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
493   71 DPR   5.132187    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
494   72  CC   5.033053    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
495   72  BH   4.651304    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
496   72 DPR   5.225125    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
497   73  CC   4.767084    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
498   73  BH   4.583107    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
499   73 DPR   5.715111    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
500   74  CC   5.125809    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
501   74  BH   4.658222    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
502   74 DPR   5.689894    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
503   75  CC   5.347350    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
504   75  BH   5.144128    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
505   75 DPR   5.774754    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
506   76  CC   5.281839    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
507   76  BH   5.167739    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
508   76 DPR   5.793248    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
509   77  CC   5.318720    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
510   77  BH   4.954516    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
511   77 DPR   6.149467    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
512   78  CC   5.361921    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
513   78  BH   5.098174    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
514   78 DPR   5.955620    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
515   79  CC   5.349554   D112A  2023-03-17    D1  12    A    7   4       5.4
516   79  CC   5.349554   D239A  2023-03-17    D2  39    A    5   7       3.9
517   79  CC   5.349554   D239B  2023-03-17    D2  39    B    4   5       5.2
518   79  CC   5.349554   F117D  2023-03-17    F1  17    D    3  13       7.2
519   79  CC   5.349554   H111B  2023-03-17    H1  11    B    7   2      12.3
520   79  CC   5.349554   H115B  2023-03-17    H1  15    B    3   3       9.7
521   79  CC   5.349554   H116A  2023-03-17    H1  16    A    4  10       9.1
522   79  CC   5.349554   H234B  2023-03-17    H2  34    B    4   7       3.4
523   79  CC   5.349554   H236C  2023-03-17    H2  36    C    2   1       6.1
524   79  CC   5.349554   J112B  2023-03-17    J1  12    B    4   4       6.3
525   79  CC   5.349554   J236C  2023-03-17    J2  36    C    5   5       7.2
526   79  CC   5.349554    L18A  2023-03-17    L1   8    A    4   6       4.2
527   79  CC   5.349554   L232A  2023-03-17    L2  32    A    5   3       5.9
528   79  CC   5.349554   L117C  2023-03-17    L1  17    C    7   3       5.4
529   79  BH   5.064459   D226B  2023-03-17    D2  26    B    5  15       5.3
530   79  BH   5.064459   D230B  2023-03-17    D2  30    B    4   3       6.1
531   79  BH   5.064459   D233B  2023-03-17    D2  33    B    6   8       5.6
532   79  BH   5.064459    D19C  2023-03-17    D1   9    C    1   8       2.6
533   79  BH   5.064459   D119D  2023-03-17    D1  19    D    5   9       2.1
534   79  BH   5.064459   D229C  2023-03-17    D2  29    C    4  13       4.3
535   79  BH   5.064459   D229D  2023-03-17    D2  29    D    2   9       4.9
536   79  BH   5.064459   D235D  2023-03-17    D2  35    D    3   6       5.3
537   79  BH   5.064459   F115A  2023-03-17    F1  15    A    6   9       8.0
538   79  BH   5.064459   F229B  2023-03-17    F2  29    B    4   5       5.9
539   79  BH   5.064459   F240B  2023-03-17    F2  40    B    5  12       5.1
540   79  BH   5.064459    F15C  2023-03-17    F1   5    C    5  26       3.0
541   79  BH   5.064459   F113C  2023-03-17    F1  13    C    3   7       4.1
542   79  BH   5.064459   F115C  2023-03-17    F1  15    C    4   4       4.5
543   79  BH   5.064459   F235C  2023-03-17    F2  35    C    5  27       5.5
544   79  BH   5.064459   F235D  2023-03-17    F2  35    D    2   1       4.6
545   79  BH   5.064459    H14B  2023-03-17    H1   4    B    2  12       5.5
546   79  BH   5.064459    H18B  2023-03-17    H1   8    B    3   9       6.5
547   79  BH   5.064459   H114A  2023-03-17    H1  14    A    5   3       6.2
548   79  BH   5.064459   H117A  2023-03-17    H1  17    A    4   6       3.6
549   79  BH   5.064459   H118B  2023-03-17    H1  18    B    5  18       5.0
550   79  BH   5.064459   H224B  2023-03-17    H2  24    B    5  29       8.5
551   79  BH   5.064459   H233B  2023-03-17    H2  33    B    4   7       3.5
552   79  BH   5.064459   H239B  2023-03-17    H2  39    B    6   2       4.0
553   79  BH   5.064459   H110C  2023-03-17    H1  10    C    6  11       5.6
554   79  BH   5.064459   H115D  2023-03-17    H1  15    D    5  13       5.6
555   79  BH   5.064459   H119C  2023-03-17    H1  19    C    7  21       6.5
556   79  BH   5.064459   H222D  2023-03-17    H2  22    D    2   3       4.4
557   79  BH   5.064459   H223C  2023-03-17    H2  23    C    7  22       4.9
558   79  BH   5.064459   H235C  2023-03-17    H2  35    C    2  13       4.8
559   79  BH   5.064459    J14B  2023-03-17    J1   4    B    2   4       4.5
560   79  BH   5.064459   J119B  2023-03-17    J1  19    B    3  11       5.3
561   79  BH   5.064459   J233A  2023-03-17    J2  33    A    3   2       4.7
562   79  BH   5.064459   J237B  2023-03-17    J2  37    B    1   7       5.3
563   79  BH   5.064459    J19C  2023-03-17    J1   9    C    5  30       8.9
564   79  BH   5.064459   J110C  2023-03-17    J1  10    C    3   1       3.9
565   79  BH   5.064459   J114D  2023-03-17    J1  14    D    4   8       3.8
566   79  BH   5.064459   J119C  2023-03-17    J1  19    C    5  20       7.1
567   79  BH   5.064459   J223D  2023-03-17    J2  23    D    4   9       5.2
568   79  BH   5.064459   J225D  2023-03-17    J2  25    D    2   5       3.0
569   79  BH   5.064459   J230C  2023-03-17    J2  30    C    4  14       8.3
570   79  BH   5.064459   J231C  2023-03-17    J2  31    C    5   6       9.2
571   79  BH   5.064459   J232D  2023-03-17    J2  32    D    3  12       4.5
572   79  BH   5.064459    L15B  2023-03-17    L1   5    B    2   6       4.2
573   79  BH   5.064459   L110A  2023-03-17    L1  10    A    2  14       2.0
574   79  BH   5.064459    L14D  2023-03-17    L1   4    D    4  10       3.8
575   79  BH   5.064459   L113C  2023-03-17    L1  13    C    3  13       3.6
576   79  BH   5.064459   L113D  2023-03-17    L1  13    D    5  22       6.5
577   79  BH   5.064459   L116C  2023-03-17    L1  16    C    6   5       7.9
578   79  BH   5.064459   L228D  2023-03-17    L2  28    D    2   7       5.1
579   79  BH   5.064459   L229D  2023-03-17    L2  29    D    3  14       3.1
580   79  BH   5.064459   L230C  2023-03-17    L2  30    C    4  11       6.6
581   79  BH   5.064459   L234C  2023-03-17    L2  34    C    5   8       4.9
582   79  BH   5.064459   L240D  2023-03-17    L2  40    D    7  13       6.0
583   79  BH   5.064459   L241D  2023-03-17    L2  41    D    4   1       4.6
584   79 DPR   6.579103   F232A  2023-03-17    F2  32    A    6   2       2.5
585   79 DPR   6.579103    F16C  2023-03-17    F1   6    C    6   9       6.6
586   79 DPR   6.579103   F240C  2023-03-17    F2  40    C    7   4      10.2
587   79 DPR   6.579103   H117B  2023-03-17    H1  17    B    5   2      17.0
588   79 DPR   6.579103    H17C  2023-03-17    H1   7    C    2  10      11.4
589   79 DPR   6.579103   H225C  2023-03-17    H2  25    C    2   6      12.4
590   79 DPR   6.579103    J16B  2023-03-17    J1   6    B    6  10       7.4
591   79 DPR   6.579103   J232B  2023-03-17    J2  32    B    6  11       4.2
592   79 DPR   6.579103   L237B  2023-03-17    L2  37    B    6   3       6.6
593   80  CC   5.606784    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
594   80  BH   5.363611    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
595   80 DPR   6.179752    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
596   81  CC   5.538654    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
597   81  BH   5.596150    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
598   81 DPR   6.430337    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
599   82  CC   5.933472    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
600   82  BH   5.403504    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
601   82 DPR   7.131757    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
602   83  CC   6.285319    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
603   83  BH   5.997065    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
604   83 DPR   7.438564    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
605   84  CC   6.162333    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
606   84  BH   5.734861    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
607   84 DPR   7.228672    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
608   85  CC   6.380242    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
609   85  BH   5.756146    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
610   85 DPR   7.325868    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
611   86  CC   6.527519   D112A  2023-03-24    D1  12    A    7   4       5.5
612   86  CC   6.527519   H111B  2023-03-24    H1  11    B    7   2      12.6
613   86  CC   6.527519   H115B  2023-03-24    H1  15    B    3   3      12.5
614   86  CC   6.527519   H116A  2023-03-24    H1  16    A    4  10      10.1
615   86  CC   6.527519   H234B  2023-03-24    H2  34    B    4   7       3.1
616   86  CC   6.527519   H236C  2023-03-24    H2  36    C    2   1       6.0
617   86  CC   6.527519   J112B  2023-03-24    J1  12    B    4   4       9.8
618   86  CC   6.527519   J236C  2023-03-24    J2  36    C    5   5       7.0
619   86  CC   6.527519    L18A  2023-03-24    L1   8    A    4   6       4.6
620   86  CC   6.527519   L232A  2023-03-24    L2  32    A    5   3       6.3
621   86  CC   6.527519   L117C  2023-03-24    L1  17    C    7   3       5.5
622   86  BH   5.915644   D226B  2023-03-24    D2  26    B    5  15       5.4
623   86  BH   5.915644   D230B  2023-03-24    D2  30    B    4   3       3.3
624   86  BH   5.915644   D233B  2023-03-24    D2  33    B    6   8       3.9
625   86  BH   5.915644   D119D  2023-03-24    D1  19    D    5   9       3.5
626   86  BH   5.915644   D229C  2023-03-24    D2  29    C    4  13       4.8
627   86  BH   5.915644   D229D  2023-03-24    D2  29    D    2   9       6.1
628   86  BH   5.915644   D235D  2023-03-24    D2  35    D    3   6       5.7
629   86  BH   5.915644   F115A  2023-03-24    F1  15    A    6   9       4.5
630   86  BH   5.915644   F229B  2023-03-24    F2  29    B    4   5       5.6
631   86  BH   5.915644   F240B  2023-03-24    F2  40    B    5  12       5.7
632   86  BH   5.915644   F115C  2023-03-24    F1  15    C    4   4       4.9
633   86  BH   5.915644   F235D  2023-03-24    F2  35    D    2   1       5.1
634   86  BH   5.915644    H14B  2023-03-24    H1   4    B    2  12       7.3
635   86  BH   5.915644    H18B  2023-03-24    H1   8    B    3   9       5.4
636   86  BH   5.915644   H114A  2023-03-24    H1  14    A    5   3       7.1
637   86  BH   5.915644   H224B  2023-03-24    H2  24    B    5  29       9.8
638   86  BH   5.915644   H233B  2023-03-24    H2  33    B    4   7       4.5
639   86  BH   5.915644   H110C  2023-03-24    H1  10    C    6  11       8.0
640   86  BH   5.915644   H115D  2023-03-24    H1  15    D    5  13       4.9
641   86  BH   5.915644   H119C  2023-03-24    H1  19    C    7  21       7.0
642   86  BH   5.915644   H222D  2023-03-24    H2  22    D    2   3       5.5
643   86  BH   5.915644   H223C  2023-03-24    H2  23    C    7  22       7.4
644   86  BH   5.915644   H235C  2023-03-24    H2  35    C    2  13       3.8
645   86  BH   5.915644    J14B  2023-03-24    J1   4    B    2   4       3.9
646   86  BH   5.915644   J119B  2023-03-24    J1  19    B    3  11       6.4
647   86  BH   5.915644   J233A  2023-03-24    J2  33    A    3   2       5.9
648   86  BH   5.915644   J237B  2023-03-24    J2  37    B    1   7       7.2
649   86  BH   5.915644    J19C  2023-03-24    J1   9    C    5  30       9.7
650   86  BH   5.915644   J110C  2023-03-24    J1  10    C    3   1       4.5
651   86  BH   5.915644   J119C  2023-03-24    J1  19    C    5  20       7.5
652   86  BH   5.915644   J223D  2023-03-24    J2  23    D    4   9       5.5
653   86  BH   5.915644   J225D  2023-03-24    J2  25    D    2   5       3.0
654   86  BH   5.915644   J230C  2023-03-24    J2  30    C    4  14       8.5
655   86  BH   5.915644   J231D  2023-03-24    J2  31    D    7  14       8.3
656   86  BH   5.915644   J232D  2023-03-24    J2  32    D    3  12       4.8
657   86  BH   5.915644    L15B  2023-03-24    L1   5    B    2   6       4.5
658   86  BH   5.915644   L110A  2023-03-24    L1  10    A    2  14       3.1
659   86  BH   5.915644    L14D  2023-03-24    L1   4    D    4  10       4.1
660   86  BH   5.915644   L113C  2023-03-24    L1  13    C    3  13       2.8
661   86  BH   5.915644   L113D  2023-03-24    L1  13    D    5  22       6.5
662   86  BH   5.915644   L228D  2023-03-24    L2  28    D    2   7       5.5
663   86  BH   5.915644   L229D  2023-03-24    L2  29    D    3  14       3.3
664   86  BH   5.915644   L230C  2023-03-24    L2  30    C    4  11       5.1
665   86  BH   5.915644   L240D  2023-03-24    L2  40    D    7  13       5.3
666   86  BH   5.915644   L241D  2023-03-24    L2  41    D    4   1       4.7
667   86 DPR   7.605860    F16C  2023-03-24    F1   6    C    6   9       9.1
668   86 DPR   7.605860   F240C  2023-03-24    F2  40    C    7   4      15.4
669   86 DPR   7.605860   H117B  2023-03-24    H1  17    B    5   2      21.7
670   86 DPR   7.605860    H17C  2023-03-24    H1   7    C    2  10      12.5
671   86 DPR   7.605860   H225C  2023-03-24    H2  25    C    2   6      18.4
672   86 DPR   7.605860    J16B  2023-03-24    J1   6    B    6  10      10.8
673   86 DPR   7.605860   J232B  2023-03-24    J2  32    B    6  11       4.4
674   86 DPR   7.605860   L237B  2023-03-24    L2  37    B    6   3       6.4
675   87  CC   6.267302    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
676   87  BH   6.294268    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
677   87 DPR   7.845873    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
678   88  CC   6.765548    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
679   88  BH   6.397516    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
680   88 DPR   8.340253    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
681   89  CC   7.029745    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
682   89  BH   6.552121    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
683   89 DPR   8.108301    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
684   90  CC   7.464191    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
685   90  BH   6.812045    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
686   90 DPR   8.834702    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
687   91  CC   7.212695    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
688   91  BH   7.195063    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
689   91 DPR   9.130418    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
690   92  CC   7.721921    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
691   92  BH   7.404194    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
692   92 DPR   9.032402    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
693   93  CC   7.563499    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
694   93  BH   7.307661    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
695   93 DPR  10.014005    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
696   94  CC   8.316324    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
697   94  BH   7.682902    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
698   94 DPR  10.302948    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
699   95  CC   8.084118    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
700   95  BH   7.805063    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
701   95 DPR  10.633471    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
702   96  CC   8.731825   H111B  2023-04-03    H1  11    B    7   2      12.1
703   96  CC   8.731825   H115B  2023-04-03    H1  15    B    3   3      14.9
704   96  CC   8.731825   H116A  2023-04-03    H1  16    A    4  10      10.2
705   96  CC   8.731825   J112B  2023-04-03    J1  12    B    4   4       7.5
706   96  CC   8.731825   J236C  2023-04-03    J2  36    C    5   5       8.3
707   96  CC   8.731825   L117C  2023-04-03    L1  17    C    7   3       6.2
708   96  BH   8.417732   D226B  2023-04-03    D2  26    B    5  15       4.9
709   96  BH   8.417732   D230B  2023-04-03    D2  30    B    4   3       3.8
710   96  BH   8.417732   D119D  2023-04-03    D1  19    D    5   9       3.2
711   96  BH   8.417732   D229C  2023-04-03    D2  29    C    4  13       4.9
712   96  BH   8.417732   D229D  2023-04-03    D2  29    D    2   9       4.1
713   96  BH   8.417732   D235D  2023-04-03    D2  35    D    3   6       6.6
714   96  BH   8.417732   F115A  2023-04-03    F1  15    A    6   9       9.9
715   96  BH   8.417732   F229B  2023-04-03    F2  29    B    4   5       5.9
716   96  BH   8.417732   F240B  2023-04-03    F2  40    B    5  12       5.6
717   96  BH   8.417732   F115C  2023-04-03    F1  15    C    4   4       5.9
718   96  BH   8.417732   F235D  2023-04-03    F2  35    D    2   1       5.5
719   96  BH   8.417732    H14B  2023-04-03    H1   4    B    2  12       0.7
720   96  BH   8.417732    H18B  2023-04-03    H1   8    B    3   9       7.1
721   96  BH   8.417732   H114A  2023-04-03    H1  14    A    5   3       8.4
722   96  BH   8.417732   H224B  2023-04-03    H2  24    B    5  29       9.4
723   96  BH   8.417732   H233B  2023-04-03    H2  33    B    4   7       4.8
724   96  BH   8.417732   H110C  2023-04-03    H1  10    C    6  11      10.6
725   96  BH   8.417732   H119C  2023-04-03    H1  19    C    7  21       6.5
726   96  BH   8.417732   H223C  2023-04-03    H2  23    C    7  22       8.4
727   96  BH   8.417732   H235C  2023-04-03    H2  35    C    2  13       3.8
728   96  BH   8.417732    J14B  2023-04-03    J1   4    B    2   4       5.1
729   96  BH   8.417732   J119B  2023-04-03    J1  19    B    3  11       5.9
730   96  BH   8.417732   J233A  2023-04-03    J2  33    A    3   2       5.8
731   96  BH   8.417732   J237B  2023-04-03    J2  37    B    1   7       4.4
732   96  BH   8.417732    J19C  2023-04-03    J1   9    C    5  30      10.2
733   96  BH   8.417732   J110C  2023-04-03    J1  10    C    3   1       5.1
734   96  BH   8.417732   J119C  2023-04-03    J1  19    C    5  20      10.1
735   96  BH   8.417732   J223D  2023-04-03    J2  23    D    4   9       5.7
736   96  BH   8.417732   J225D  2023-04-03    J2  25    D    2   5       6.3
737   96  BH   8.417732   J230C  2023-04-03    J2  30    C    4  14       9.4
738   96  BH   8.417732   J231D  2023-04-03    J2  31    D    7  14      12.1
739   96  BH   8.417732   J232D  2023-04-03    J2  32    D    3  12       5.2
740   96  BH   8.417732    L15B  2023-04-03    L1   5    B    2   6       3.9
741   96  BH   8.417732   L110A  2023-04-03    L1  10    A    2  14       3.6
742   96  BH   8.417732    L14D  2023-04-03    L1   4    D    4  10       7.1
743   96  BH   8.417732   L113C  2023-04-03    L1  13    C    3  13       3.6
744   96  BH   8.417732   L113D  2023-04-03    L1  13    D    5  22       7.3
745   96  BH   8.417732   L228D  2023-04-03    L2  28    D    2   7       7.1
746   96  BH   8.417732   L229D  2023-04-03    L2  29    D    3  14       4.1
747   96  BH   8.417732   L230C  2023-04-03    L2  30    C    4  11       4.3
748   96  BH   8.417732   L240D  2023-04-03    L2  40    D    7  13       6.6
749   96  BH   8.417732   L241D  2023-04-03    L2  41    D    4   1       5.8
750   96 DPR  11.182425    F16C  2023-04-03    F1   6    C    6   9      15.4
751   96 DPR  11.182425   F240C  2023-04-03    F2  40    C    7   4      22.5
752   96 DPR  11.182425   H117B  2023-04-03    H1  17    B    5   2      25.9
753   96 DPR  11.182425    H17C  2023-04-03    H1   7    C    2  10      13.3
754   96 DPR  11.182425   H225C  2023-04-03    H2  25    C    2   6      27.7
755   97  CC   8.836426    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
756   97  BH   8.288714    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
757   97 DPR  11.757594    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
758   98  CC   9.267155    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
759   98  BH   8.162741    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
760   98 DPR  11.918344    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
761   99  CC   9.435969    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
762   99  BH   8.967746    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
763   99 DPR  12.570296    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
764  100  CC   9.755982    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
765  100  BH   8.901483    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
766  100 DPR  13.015052    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
767  101  CC  10.426101    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
768  101  BH   9.321250    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
769  101 DPR  13.644455    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
770  102  CC  10.538790    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
771  102  BH   9.562548    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
772  102 DPR  13.952632    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
773  103  CC  11.008650   H111B  2023-04-10    H1  11    B    7   2      16.2
774  103  CC  11.008650   H115B  2023-04-10    H1  15    B    3   3      23.8
775  103  CC  11.008650   H116A  2023-04-10    H1  16    A    4  10      12.3
776  103  CC  11.008650   J112B  2023-04-10    J1  12    B    4   4      10.1
777  103  CC  11.008650   J236C  2023-04-10    J2  36    C    5   5       5.9
778  103  BH  10.021828   D226B  2023-04-10    D2  26    B    5  15       5.7
779  103  BH  10.021828   D230B  2023-04-10    D2  30    B    4   3       6.8
780  103  BH  10.021828   D119D  2023-04-10    D1  19    D    5   9       3.6
781  103  BH  10.021828   D229C  2023-04-10    D2  29    C    4  13       5.8
782  103  BH  10.021828   D229D  2023-04-10    D2  29    D    2   9       4.3
783  103  BH  10.021828   D235D  2023-04-10    D2  35    D    3   6       8.4
784  103  BH  10.021828   F115A  2023-04-10    F1  15    A    6   9      17.4
785  103  BH  10.021828   F229B  2023-04-10    F2  29    B    4   5       7.2
786  103  BH  10.021828   F240B  2023-04-10    F2  40    B    5  12       6.6
787  103  BH  10.021828   F115C  2023-04-10    F1  15    C    4   4       5.5
788  103  BH  10.021828   F235D  2023-04-10    F2  35    D    2   1       5.7
789  103  BH  10.021828    H14B  2023-04-10    H1   4    B    2  12       6.1
790  103  BH  10.021828    H18B  2023-04-10    H1   8    B    3   9       8.8
791  103  BH  10.021828   H114A  2023-04-10    H1  14    A    5   3      12.3
792  103  BH  10.021828   H224B  2023-04-10    H2  24    B    5  29      13.4
793  103  BH  10.021828   H233B  2023-04-10    H2  33    B    4   7       5.7
794  103  BH  10.021828   H110C  2023-04-10    H1  10    C    6  11      13.2
795  103  BH  10.021828   H119C  2023-04-10    H1  19    C    7  21       5.1
796  103  BH  10.021828   H223C  2023-04-10    H2  23    C    7  22       8.6
797  103  BH  10.021828   H235C  2023-04-10    H2  35    C    2  13       3.9
798  103  BH  10.021828    J14B  2023-04-10    J1   4    B    2   4       6.0
799  103  BH  10.021828   J119B  2023-04-10    J1  19    B    3  11       6.9
800  103  BH  10.021828   J233A  2023-04-10    J2  33    A    3   2       6.4
801  103  BH  10.021828   J237B  2023-04-10    J2  37    B    1   7      12.2
802  103  BH  10.021828    J19C  2023-04-10    J1   9    C    5  30       8.5
803  103  BH  10.021828   J110C  2023-04-10    J1  10    C    3   1       7.0
804  103  BH  10.021828   J119C  2023-04-10    J1  19    C    5  20      11.4
805  103  BH  10.021828   J225D  2023-04-10    J2  25    D    2   5       7.7
806  103  BH  10.021828   J230C  2023-04-10    J2  30    C    4  14       9.0
807  103  BH  10.021828   J231D  2023-04-10    J2  31    D    7  14      10.2
808  103  BH  10.021828   J232D  2023-04-10    J2  32    D    3  12       6.0
809  103  BH  10.021828   L110A  2023-04-10    L1  10    A    2  14       2.1
810  103  BH  10.021828    L14D  2023-04-10    L1   4    D    4  10       4.1
811  103  BH  10.021828   L113C  2023-04-10    L1  13    C    3  13       2.5
812  103  BH  10.021828   L113D  2023-04-10    L1  13    D    5  22       9.0
813  103  BH  10.021828   L228D  2023-04-10    L2  28    D    2   7       6.3
814  103  BH  10.021828   L229D  2023-04-10    L2  29    D    3  14       3.6
815  103  BH  10.021828   L240D  2023-04-10    L2  40    D    7  13       8.4
816  103  BH  10.021828   L241D  2023-04-10    L2  41    D    4   1       4.4
817  103 DPR  14.575802    F16C  2023-04-10    F1   6    C    6   9      16.6
818  103 DPR  14.575802   H117B  2023-04-10    H1  17    B    5   2      38.2
819  103 DPR  14.575802   H225C  2023-04-10    H2  25    C    2   6      37.4
820  103 DPR  14.575802    J16B  2023-04-10    J1   6    B    6  10      23.6
821  104  CC  11.508182    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
822  104  BH  10.489636    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
823  104 DPR  15.313637    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
824  105  CC  11.898376    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
825  105  BH  10.626219    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
826  105 DPR  15.851956    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
827  106  CC  12.479055    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
828  106  BH  11.074797    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
829  106 DPR  16.422490    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
830  107  CC  12.644022    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
831  107  BH  11.206877    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
832  107 DPR  16.782215    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
833  108  CC  13.268446    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
834  108  BH  11.353314    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
835  108 DPR  17.731581    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
836  109  CC  13.168846    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
837  109  BH  12.387764    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
838  109 DPR  18.272415    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
839  110  CC  13.816665   H111B  2023-04-17    H1  11    B    7   2      14.4
840  110  CC  13.816665   H115B  2023-04-17    H1  15    B    3   3      32.6
841  110  CC  13.816665   H116A  2023-04-17    H1  16    A    4  10      15.1
842  110  CC  13.816665   J112B  2023-04-17    J1  12    B    4   4       9.0
843  110  BH  12.390678   D226B  2023-04-17    D2  26    B    5  15       7.9
844  110  BH  12.390678   D230B  2023-04-17    D2  30    B    4   3       8.1
845  110  BH  12.390678   D119D  2023-04-17    D1  19    D    5   9       3.3
846  110  BH  12.390678   D229C  2023-04-17    D2  29    C    4  13       6.8
847  110  BH  12.390678   D229D  2023-04-17    D2  29    D    2   9       7.1
848  110  BH  12.390678   D235D  2023-04-17    D2  35    D    3   6      10.2
849  110  BH  12.390678   F115A  2023-04-17    F1  15    A    6   9      22.9
850  110  BH  12.390678   F229B  2023-04-17    F2  29    B    4   5       8.2
851  110  BH  12.390678   F240B  2023-04-17    F2  40    B    5  12       8.0
852  110  BH  12.390678   F235D  2023-04-17    F2  35    D    2   1       7.3
853  110  BH  12.390678    H14B  2023-04-17    H1   4    B    2  12      11.4
854  110  BH  12.390678    H18B  2023-04-17    H1   8    B    3   9      15.2
855  110  BH  12.390678   H114A  2023-04-17    H1  14    A    5   3      13.1
856  110  BH  12.390678   H224B  2023-04-17    H2  24    B    5  29      19.2
857  110  BH  12.390678   H233B  2023-04-17    H2  33    B    4   7       6.2
858  110  BH  12.390678   H110C  2023-04-17    H1  10    C    6  11      18.3
859  110  BH  12.390678   H223C  2023-04-17    H2  23    C    7  22      11.8
860  110  BH  12.390678    J14B  2023-04-17    J1   4    B    2   4      12.9
861  110  BH  12.390678   J119B  2023-04-17    J1  19    B    3  11      12.0
862  110  BH  12.390678   J233A  2023-04-17    J2  33    A    3   2       9.0
863  110  BH  12.390678   J237B  2023-04-17    J2  37    B    1   7      18.4
864  110  BH  12.390678    J19C  2023-04-17    J1   9    C    5  30       9.8
865  110  BH  12.390678   J110C  2023-04-17    J1  10    C    3   1      13.4
866  110  BH  12.390678   J119C  2023-04-17    J1  19    C    5  20      18.4
867  110  BH  12.390678   J223D  2023-04-17    J2  23    D    4   9       6.1
868  110  BH  12.390678   J225D  2023-04-17    J2  25    D    2   5       8.8
869  110  BH  12.390678   J230C  2023-04-17    J2  30    C    4  14      13.1
870  110  BH  12.390678   J231D  2023-04-17    J2  31    D    7  14      11.2
871  110  BH  12.390678   J232D  2023-04-17    J2  32    D    3  12      10.8
872  110  BH  12.390678   L110A  2023-04-17    L1  10    A    2  14      11.1
873  110  BH  12.390678    L14D  2023-04-17    L1   4    D    4  10      14.4
874  110  BH  12.390678   L113C  2023-04-17    L1  13    C    3  13      16.0
875  110  BH  12.390678   L113D  2023-04-17    L1  13    D    5  22      12.4
876  110  BH  12.390678   L228D  2023-04-17    L2  28    D    2   7       7.1
877  110  BH  12.390678   L229D  2023-04-17    L2  29    D    3  14       5.6
878  110  BH  12.390678   L240D  2023-04-17    L2  40    D    7  13      11.8
879  110  BH  12.390678   L241D  2023-04-17    L2  41    D    4   1       6.9
880  110 DPR  19.242244    F16C  2023-04-17    F1   6    C    6   9      16.1
881  110 DPR  19.242244   H117B  2023-04-17    H1  17    B    5   2      42.7
882  110 DPR  19.242244   H225C  2023-04-17    H2  25    C    2   6      36.2
883  110 DPR  19.242244    J16B  2023-04-17    J1   6    B    6  10      24.5
884  111  CC  14.880424    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
885  111  BH  12.881263    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
886  111 DPR  19.837310    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
887  112  CC  15.068294    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
888  112  BH  13.266819    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
889  112 DPR  21.195502    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
890  113  CC  15.417283    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
891  113  BH  13.784213    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
892  113 DPR  21.746087    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
893  114  CC  16.159700    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
894  114  BH  14.229887    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
895  114 DPR  22.785871    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
896  115  CC  16.944652    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
897  115  BH  14.995171    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
898  115 DPR  23.644658    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
899  116  CC  17.345762    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
900  116  BH  15.547021    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
901  116 DPR  24.884839    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
902  117  CC  17.804378   H111B  2023-04-24    H1  11    B    7   2      27.1
903  117  CC  17.804378   H115B  2023-04-24    H1  15    B    3   3      37.0
904  117  CC  17.804378   H116A  2023-04-24    H1  16    A    4  10      15.9
905  117  BH  15.891774   D226B  2023-04-24    D2  26    B    5  15      15.9
906  117  BH  15.891774   D230B  2023-04-24    D2  30    B    4   3      12.6
907  117  BH  15.891774   D119D  2023-04-24    D1  19    D    5   9       2.1
908  117  BH  15.891774   D229C  2023-04-24    D2  29    C    4  13       8.1
909  117  BH  15.891774   D229D  2023-04-24    D2  29    D    2   9      10.9
910  117  BH  15.891774   D235D  2023-04-24    D2  35    D    3   6      16.3
911  117  BH  15.891774   F115A  2023-04-24    F1  15    A    6   9      33.1
912  117  BH  15.891774   F229B  2023-04-24    F2  29    B    4   5       9.1
913  117  BH  15.891774   F240B  2023-04-24    F2  40    B    5  12      10.1
914  117  BH  15.891774   F235D  2023-04-24    F2  35    D    2   1       1.8
915  117  BH  15.891774    H14B  2023-04-24    H1   4    B    2  12      14.9
916  117  BH  15.891774    H18B  2023-04-24    H1   8    B    3   9      21.4
917  117  BH  15.891774   H114A  2023-04-24    H1  14    A    5   3      17.5
918  117  BH  15.891774   H224B  2023-04-24    H2  24    B    5  29      28.3
919  117  BH  15.891774   H233B  2023-04-24    H2  33    B    4   7      13.7
920  117  BH  15.891774   H110C  2023-04-24    H1  10    C    6  11      24.4
921  117  BH  15.891774   H223C  2023-04-24    H2  23    C    7  22      12.2
922  117  BH  15.891774    J14B  2023-04-24    J1   4    B    2   4       5.2
923  117  BH  15.891774   J119B  2023-04-24    J1  19    B    3  11      11.9
924  117  BH  15.891774   J233A  2023-04-24    J2  33    A    3   2       7.9
925  117  BH  15.891774   J237B  2023-04-24    J2  37    B    1   7      25.3
926  117  BH  15.891774   J110C  2023-04-24    J1  10    C    3   1      14.8
927  117  BH  15.891774   J119C  2023-04-24    J1  19    C    5  20      28.3
928  117  BH  15.891774   J223D  2023-04-24    J2  23    D    4   9       8.4
929  117  BH  15.891774   J225D  2023-04-24    J2  25    D    2   5      14.1
930  117  BH  15.891774   J230C  2023-04-24    J2  30    C    4  14      15.3
931  117  BH  15.891774   J231D  2023-04-24    J2  31    D    7  14      10.6
932  117  BH  15.891774   J232D  2023-04-24    J2  32    D    3  12      19.2
933  117  BH  15.891774   L110A  2023-04-24    L1  10    A    2  14      14.4
934  117  BH  15.891774    L14D  2023-04-24    L1   4    D    4  10      20.6
935  117  BH  15.891774   L113C  2023-04-24    L1  13    C    3  13      27.7
936  117  BH  15.891774   L113D  2023-04-24    L1  13    D    5  22      17.5
937  117  BH  15.891774   L228D  2023-04-24    L2  28    D    2   7       6.9
938  117  BH  15.891774   L229D  2023-04-24    L2  29    D    3  14      11.8
939  117  BH  15.891774   L240D  2023-04-24    L2  40    D    7  13      14.1
940  117  BH  15.891774   L241D  2023-04-24    L2  41    D    4   1      12.3
941  117 DPR  25.382756    F16C  2023-04-24    F1   6    C    6   9      16.0
942  117 DPR  25.382756   H117B  2023-04-24    H1  17    B    5   2      42.3
943  117 DPR  25.382756   H225C  2023-04-24    H2  25    C    2   6      36.1
944  117 DPR  25.382756    J16B  2023-04-24    J1   6    B    6  10      26.7
945  118  CC  18.755423    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
946  118  BH  16.526043    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
947  118 DPR  26.452630    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
948  119  CC  18.918602    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
949  119  BH  17.068545    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
950  119 DPR  27.590567    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
951  120  CC  20.082350    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
952  120  BH  17.496580    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
953  120 DPR  28.496450    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
954  121  CC  20.496085    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
955  121  BH  18.333086    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
956  121 DPR  29.730818    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
957  122  CC  21.371723    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
958  122  BH  19.145254    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
959  122 DPR  30.742783    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
960  123  CC  22.219562    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
961  123  BH  19.626691    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
962  123 DPR  32.191712    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
963  124  CC  23.154456   H111B  2023-05-01    H1  11    B    7   2      15.5
964  124  CC  23.154456   H115B  2023-05-01    H1  15    B    3   3      34.7
965  124  CC  23.154456   H116A  2023-05-01    H1  16    A    4  10      13.5
966  124  BH  20.010639   D226B  2023-05-01    D2  26    B    5  15      20.5
967  124  BH  20.010639   D230B  2023-05-01    D2  30    B    4   3      17.2
968  124  BH  20.010639   D119D  2023-05-01    D1  19    D    5   9       3.1
969  124  BH  20.010639   D229C  2023-05-01    D2  29    C    4  13      10.8
970  124  BH  20.010639   D229D  2023-05-01    D2  29    D    2   9      13.1
971  124  BH  20.010639   D235D  2023-05-01    D2  35    D    3   6       3.9
972  124  BH  20.010639   F115A  2023-05-01    F1  15    A    6   9      42.5
973  124  BH  20.010639   F229B  2023-05-01    F2  29    B    4   5      12.5
974  124  BH  20.010639   F240B  2023-05-01    F2  40    B    5  12      16.5
975  124  BH  20.010639   F235D  2023-05-01    F2  35    D    2   1      10.2
976  124  BH  20.010639    H14B  2023-05-01    H1   4    B    2  12      18.6
977  124  BH  20.010639    H18B  2023-05-01    H1   8    B    3   9      30.5
978  124  BH  20.010639   H114A  2023-05-01    H1  14    A    5   3      19.5
979  124  BH  20.010639   H224B  2023-05-01    H2  24    B    5  29      38.4
980  124  BH  20.010639   H233B  2023-05-01    H2  33    B    4   7      19.1
981  124  BH  20.010639   H110C  2023-05-01    H1  10    C    6  11      29.9
982  124  BH  20.010639   H223C  2023-05-01    H2  23    C    7  22      12.1
983  124  BH  20.010639    J14B  2023-05-01    J1   4    B    2   4      11.5
984  124  BH  20.010639   J119B  2023-05-01    J1  19    B    3  11      18.4
985  124  BH  20.010639   J233A  2023-05-01    J2  33    A    3   2       9.5
986  124  BH  20.010639   J237B  2023-05-01    J2  37    B    1   7      33.0
987  124  BH  20.010639   J110C  2023-05-01    J1  10    C    3   1      31.5
988  124  BH  20.010639   J119C  2023-05-01    J1  19    C    5  20      44.3
989  124  BH  20.010639   J223D  2023-05-01    J2  23    D    4   9      10.2
990  124  BH  20.010639   J225D  2023-05-01    J2  25    D    2   5      17.0
991  124  BH  20.010639   J230C  2023-05-01    J2  30    C    4  14      23.5
992  124  BH  20.010639   J231D  2023-05-01    J2  31    D    7  14      10.5
993  124  BH  20.010639   J232D  2023-05-01    J2  32    D    3  12      23.5
994  124  BH  20.010639   L110A  2023-05-01    L1  10    A    2  14      21.0
995  124  BH  20.010639    L14D  2023-05-01    L1   4    D    4  10      24.7
996  124  BH  20.010639   L113C  2023-05-01    L1  13    C    3  13      50.6
997  124  BH  20.010639   L113D  2023-05-01    L1  13    D    5  22      21.1
998  124  BH  20.010639   L228D  2023-05-01    L2  28    D    2   7       7.5
999  124  BH  20.010639   L229D  2023-05-01    L2  29    D    3  14      16.5
1000 124  BH  20.010639   L240D  2023-05-01    L2  40    D    7  13      15.0
1001 124  BH  20.010639   L241D  2023-05-01    L2  41    D    4   1      16.5
1002 124 DPR  33.187743    F16C  2023-05-01    F1   6    C    6   9      17.0
1003 124 DPR  33.187743   H117B  2023-05-01    H1  17    B    5   2      42.4
1004 124 DPR  33.187743    J16B  2023-05-01    J1   6    B    6  10      23.5
1005 125  CC  24.043395    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1006 125  BH  20.739748    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1007 125 DPR  34.506160    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1008 126  CC  24.714385    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1009 126  BH  21.679890    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1010 126 DPR  35.512855    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1011 127  CC  25.573105    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1012 127  BH  22.111137    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1013 127 DPR  36.744408    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1014 128  CC  26.000498    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1015 128  BH  23.349049    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1016 128 DPR  37.941599    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1017 129  CC  27.037045    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1018 129  BH  23.860166    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1019 129 DPR  39.503336    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1020 130  CC  27.967606    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1021 130  BH  25.097028    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1022 130 DPR  40.693179    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1023 131  CC  28.612470   H116A  2023-05-08    H1  16    A    4  10      15.1
1024 131  BH  25.479420   D226B  2023-05-08    D2  26    B    5  15      23.1
1025 131  BH  25.479420   D230B  2023-05-08    D2  30    B    4   3      22.1
1026 131  BH  25.479420   D119D  2023-05-08    D1  19    D    5   9       3.1
1027 131  BH  25.479420   D229C  2023-05-08    D2  29    C    4  13      13.3
1028 131  BH  25.479420   D229D  2023-05-08    D2  29    D    2   9      15.8
1029 131  BH  25.479420   D235D  2023-05-08    D2  35    D    3   6      35.2
1030 131  BH  25.479420   F115A  2023-05-08    F1  15    A    6   9      51.3
1031 131  BH  25.479420   F229B  2023-05-08    F2  29    B    4   5      14.6
1032 131  BH  25.479420   F240B  2023-05-08    F2  40    B    5  12      18.7
1033 131  BH  25.479420   F235D  2023-05-08    F2  35    D    2   1      12.2
1034 131  BH  25.479420    H14B  2023-05-08    H1   4    B    2  12      20.8
1035 131  BH  25.479420    H18B  2023-05-08    H1   8    B    3   9      36.1
1036 131  BH  25.479420   H114A  2023-05-08    H1  14    A    5   3      21.6
1037 131  BH  25.479420   H224B  2023-05-08    H2  24    B    5  29      48.3
1038 131  BH  25.479420   H233B  2023-05-08    H2  33    B    4   7      20.8
1039 131  BH  25.479420   H110C  2023-05-08    H1  10    C    6  11      34.0
1040 131  BH  25.479420   H223C  2023-05-08    H2  23    C    7  22      12.1
1041 131  BH  25.479420    J14B  2023-05-08    J1   4    B    2   4       4.4
1042 131  BH  25.479420   J119B  2023-05-08    J1  19    B    3  11      18.6
1043 131  BH  25.479420   J233A  2023-05-08    J2  33    A    3   2       5.7
1044 131  BH  25.479420   J237B  2023-05-08    J2  37    B    1   7      39.2
1045 131  BH  25.479420   J110C  2023-05-08    J1  10    C    3   1      33.1
1046 131  BH  25.479420   J119C  2023-05-08    J1  19    C    5  20      60.5
1047 131  BH  25.479420   J223D  2023-05-08    J2  23    D    4   9      13.2
1048 131  BH  25.479420   J225D  2023-05-08    J2  25    D    2   5      20.4
1049 131  BH  25.479420   J230C  2023-05-08    J2  30    C    4  14      24.6
1050 131  BH  25.479420   J231D  2023-05-08    J2  31    D    7  14       8.9
1051 131  BH  25.479420   J232D  2023-05-08    J2  32    D    3  12      24.0
1052 131  BH  25.479420   L110A  2023-05-08    L1  10    A    2  14      23.3
1053 131  BH  25.479420    L14D  2023-05-08    L1   4    D    4  10      27.3
1054 131  BH  25.479420   L113C  2023-05-08    L1  13    C    3  13      65.9
1055 131  BH  25.479420   L113D  2023-05-08    L1  13    D    5  22      21.2
1056 131  BH  25.479420   L228D  2023-05-08    L2  28    D    2   7       5.0
1057 131  BH  25.479420   L229D  2023-05-08    L2  29    D    3  14      17.8
1058 131  BH  25.479420   L240D  2023-05-08    L2  40    D    7  13      14.1
1059 131  BH  25.479420   L241D  2023-05-08    L2  41    D    4   1      11.5
1060 131 DPR  42.773905   H117B  2023-05-08    H1  17    B    5   2      42.2
1061 132  CC  30.088978    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1062 132  BH  26.186182    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1063 132 DPR  43.647173    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1064 133  CC  31.196331    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1065 133  BH  26.816032    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1066 133 DPR  45.194823    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1067 134  CC  31.765057    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1068 134  BH  27.684059    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1069 134 DPR  47.016312    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1070 135  CC  32.769783    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1071 135  BH  28.603156    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1072 135 DPR  48.455928    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1073 136  CC  34.377848    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1074 136  BH  29.351422    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1075 136 DPR  49.585429    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1076 137  CC  34.846974    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1077 137  BH  30.192838    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1078 137 DPR  51.356972    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1079 138  CC  36.207839    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1080 138  BH  31.171344    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1081 138 DPR  52.890046    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1082 139  CC  36.894982    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1083 139  BH  31.797720    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1084 139 DPR  54.586425    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1085 140  CC  38.179100    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1086 140  BH  33.129145    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1087 140 DPR  56.109876    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1088 141  CC  39.179105    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1089 141  BH  33.856948    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1090 141 DPR  57.547399    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1091 142  CC  40.066411    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1092 142  BH  34.647191    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1093 142 DPR  59.137869    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1094 143  CC  41.408780    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1095 143  BH  35.969600    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1096 143 DPR  61.069291    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1097 144  CC  42.410637    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1098 144  BH  36.518004    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1099 144 DPR  62.757482    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1100 145  CC  43.255354    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1101 145  BH  37.514377    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1102 145 DPR  63.820550    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1103 146  CC  44.361852    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1104 146  BH  38.544644    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1105 146 DPR  65.909251    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1106 147  CC  45.766500    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1107 147  BH  39.135822    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1108 147 DPR  67.240561    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1109 148  CC  46.641047    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1110 148  BH  40.452533    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1111 148 DPR  69.318113    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1112 149  CC  47.995462    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1113 149  BH  41.020525    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1114 149 DPR  70.643228    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1115 150  CC  48.851410    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1116 150  BH  41.902647    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1117 150 DPR  72.283466    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1118 151  CC  49.380357    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1119 151  BH  42.724570    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1120 151 DPR  74.147447    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1121 152  CC  50.739702    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1122 152  BH  44.022608    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1123 152 DPR  75.446032    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1124 153  CC  52.095573    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1125 153  BH  44.791203    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1126 153 DPR  76.931500    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1127 154  CC  53.059480    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1128 154  BH  45.699247    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1129 154 DPR  78.687504    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1130 155  CC  54.122118    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1131 155  BH  46.616206    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1132 155 DPR  80.313658    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1133 156  CC  55.176722    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1134 156  BH  47.448632    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1135 156 DPR  81.731053    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1136 157  CC  56.206381    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1137 157  BH  48.023528    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1138 157 DPR  83.232497    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1139 158  CC  57.065209    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1140 158  BH  49.060973    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1141 158 DPR  84.537921    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1142 159  CC  57.967421    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1143 159  BH  49.686407   D226B  2023-06-05    D2  26    B    5  15      77.4
1144 159  BH  49.686407   D230B  2023-06-05    D2  30    B    4   3      52.1
1145 159  BH  49.686407   D229C  2023-06-05    D2  29    C    4  13      23.7
1146 159  BH  49.686407   D229D  2023-06-05    D2  29    D    2   9      34.5
1147 159  BH  49.686407   D235D  2023-06-05    D2  35    D    3   6      80.5
1148 159  BH  49.686407   F115A  2023-06-05    F1  15    A    6   9      88.9
1149 159  BH  49.686407   F229B  2023-06-05    F2  29    B    4   5      26.3
1150 159  BH  49.686407   F240B  2023-06-05    F2  40    B    5  12      47.0
1151 159  BH  49.686407   F235D  2023-06-05    F2  35    D    2   1      37.5
1152 159  BH  49.686407    H14B  2023-06-05    H1   4    B    2  12      28.2
1153 159  BH  49.686407    H18B  2023-06-05    H1   8    B    3   9      76.1
1154 159  BH  49.686407   H114A  2023-06-05    H1  14    A    5   3      19.9
1155 159  BH  49.686407   H224B  2023-06-05    H2  24    B    5  29      70.8
1156 159  BH  49.686407   H233B  2023-06-05    H2  33    B    4   7      21.1
1157 159  BH  49.686407   H110C  2023-06-05    H1  10    C    6  11      44.8
1158 159  BH  49.686407   J119B  2023-06-05    J1  19    B    3  11      31.9
1159 159  BH  49.686407   J237B  2023-06-05    J2  37    B    1   7      39.9
1160 159  BH  49.686407   J110C  2023-06-05    J1  10    C    3   1      70.0
1161 159  BH  49.686407   J119C  2023-06-05    J1  19    C    5  20     102.4
1162 159  BH  49.686407   J223D  2023-06-05    J2  23    D    4   9      51.4
1163 159  BH  49.686407   J225D  2023-06-05    J2  25    D    2   5      71.0
1164 159  BH  49.686407   J230C  2023-06-05    J2  30    C    4  14      32.3
1165 159  BH  49.686407   J232D  2023-06-05    J2  32    D    3  12      27.8
1166 159  BH  49.686407   L110A  2023-06-05    L1  10    A    2  14      61.7
1167 159  BH  49.686407    L14D  2023-06-05    L1   4    D    4  10      77.6
1168 159  BH  49.686407   L113C  2023-06-05    L1  13    C    3  13     120.6
1169 159  BH  49.686407   L113D  2023-06-05    L1  13    D    5  22      58.9
1170 159  BH  49.686407   L229D  2023-06-05    L2  29    D    3  14      31.4
1171 159  BH  49.686407   L240D  2023-06-05    L2  40    D    7  13      19.7
1172 159  BH  49.686407   L241D  2023-06-05    L2  41    D    4   1      49.4
1173 159 DPR  86.186067    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1174 160  CC  59.041226    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1175 160  BH  50.350054    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1176 160 DPR  87.555664    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1177 161  CC  59.742334    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1178 161  BH  51.367100    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1179 161 DPR  89.345603    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1180 162  CC  60.512758    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1181 162  BH  52.269627    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1182 162 DPR  90.311694    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1183 163  CC  61.477930    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1184 163  BH  53.101869    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1185 163 DPR  91.956132    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1186 164  CC  62.601635    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1187 164  BH  53.600240    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1188 164 DPR  93.283346    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1189 165  CC  63.096784    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1190 165  BH  53.955111    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1191 165 DPR  94.238881    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1192 166  CC  63.860656    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1193 166  BH  55.152038    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1194 166 DPR  95.663059    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1195 167  CC  64.749152    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1196 167  BH  55.594440    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1197 167 DPR  96.819425    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1198 168  CC  65.679603    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1199 168  BH  56.194431    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1200 168 DPR  98.041667    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1201 169  CC  65.924640    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1202 169  BH  56.886229    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1203 169 DPR  99.156191    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1204 170  CC  66.847329    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1205 170  BH  57.547756    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1206 170 DPR  99.976718    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1207 171  CC  67.554956    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1208 171  BH  57.826521    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1209 171 DPR 101.064611    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1210 172  CC  68.358637    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1211 172  BH  58.454705    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1212 172 DPR 101.892991    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1213 173  CC  69.080368    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1214 173  BH  58.973871    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1215 173 DPR 102.847723    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1216 174  CC  69.430609    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1217 174  BH  59.606107    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1218 174 DPR 103.688635    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1219 175  CC  69.789069    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1220 175  BH  59.810451    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1221 175 DPR 104.570287    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1222 176  CC  70.640889    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1223 176  BH  60.507520    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1224 176 DPR 105.639181    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1225 177  CC  70.960184    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1226 177  BH  60.571807    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1227 177 DPR 106.039139    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1228 178  CC  71.352954    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1229 178  BH  61.017703    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1230 178 DPR 106.692398    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1231 179  CC  71.624029    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1232 179  BH  61.353802    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1233 179 DPR 107.253589    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1234 180  CC  72.217475    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1235 180  BH  61.772673    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1236 180 DPR 107.785244    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1237 181  CC  72.830166    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1238 181  BH  62.252335    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1239 181 DPR 108.556509    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1240 182  CC  72.873571    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1241 182  BH  62.449706    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1242 182 DPR 109.040189    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1243 183  CC  73.027912    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1244 183  BH  62.638459    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1245 183 DPR 109.570341    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1246 184  CC  73.385856    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1247 184  BH  62.537297    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1248 184 DPR 110.053257    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1249 185  CC  74.004917    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1250 185  BH  63.253726    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1251 185 DPR 110.318259    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1252 186  CC  74.164970    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1253 186  BH  63.523171    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1254 186 DPR 110.604714    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1255 187  CC  74.629540    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1256 187  BH  63.677912   D226B  2023-07-03    D2  26    B    5  15      67.4
1257 187  BH  63.677912   D230B  2023-07-03    D2  30    B    4   3      76.3
1258 187  BH  63.677912   D229C  2023-07-03    D2  29    C    4  13      26.9
1259 187  BH  63.677912   D229D  2023-07-03    D2  29    D    2   9      42.7
1260 187  BH  63.677912   D235D  2023-07-03    D2  35    D    3   6      83.2
1261 187  BH  63.677912   F115A  2023-07-03    F1  15    A    6   9     110.4
1262 187  BH  63.677912   F229B  2023-07-03    F2  29    B    4   5      34.3
1263 187  BH  63.677912   F240B  2023-07-03    F2  40    B    5  12      70.9
1264 187  BH  63.677912   F235D  2023-07-03    F2  35    D    2   1      67.8
1265 187  BH  63.677912    H18B  2023-07-03    H1   8    B    3   9      81.6
1266 187  BH  63.677912   H114A  2023-07-03    H1  14    A    5   3      18.7
1267 187  BH  63.677912   H224B  2023-07-03    H2  24    B    5  29      69.0
1268 187  BH  63.677912   H233B  2023-07-03    H2  33    B    4   7      18.0
1269 187  BH  63.677912   H110C  2023-07-03    H1  10    C    6  11      45.1
1270 187  BH  63.677912   J119B  2023-07-03    J1  19    B    3  11      32.1
1271 187  BH  63.677912   J110C  2023-07-03    J1  10    C    3   1      67.0
1272 187  BH  63.677912   J119C  2023-07-03    J1  19    C    5  20     101.0
1273 187  BH  63.677912   J223D  2023-07-03    J2  23    D    4   9      49.0
1274 187  BH  63.677912   J225D  2023-07-03    J2  25    D    2   5      82.0
1275 187  BH  63.677912   J230C  2023-07-03    J2  30    C    4  14      35.0
1276 187  BH  63.677912   J232D  2023-07-03    J2  32    D    3  12      29.0
1277 187  BH  63.677912   L110A  2023-07-03    L1  10    A    2  14      89.5
1278 187  BH  63.677912    L14D  2023-07-03    L1   4    D    4  10      92.5
1279 187  BH  63.677912   L113C  2023-07-03    L1  13    C    3  13      82.0
1280 187  BH  63.677912   L113D  2023-07-03    L1  13    D    5  22     113.0
1281 187  BH  63.677912   L229D  2023-07-03    L2  29    D    3  14      28.0
1282 187  BH  63.677912   L240D  2023-07-03    L2  40    D    7  13      17.0
1283 187  BH  63.677912   L241D  2023-07-03    L2  41    D    4   1      67.5
1284 187 DPR 111.193817    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1285 188  CC  74.716680    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1286 188  BH  63.795310    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1287 188 DPR 111.588816    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1288 189  CC  74.964410    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1289 189  BH  63.961600    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1290 189 DPR 111.543074    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1291 190  CC  74.934354    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1292 190  BH  64.305275    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1293 190 DPR 112.154336    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1294 191  CC  75.113608    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1295 191  BH  64.426350    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1296 191 DPR 112.336997    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1297 192  CC  75.304927    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1298 192  BH  64.149685    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1299 192 DPR 112.480630    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1300 193  CC  75.477320    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1301 193  BH  64.410890    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1302 193 DPR 112.485540    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1303 194  CC  75.502619    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1304 194  BH  64.589675   D226B  2023-07-10    D2  26    B    5  15      70.0
1305 194  BH  64.589675   D230B  2023-07-10    D2  30    B    4   3      73.0
1306 194  BH  64.589675   D229C  2023-07-10    D2  29    C    4  13      26.5
1307 194  BH  64.589675   D229D  2023-07-10    D2  29    D    2   9      43.0
1308 194  BH  64.589675   D235D  2023-07-10    D2  35    D    3   6      85.8
1309 194  BH  64.589675   F115A  2023-07-10    F1  15    A    6   9     113.2
1310 194  BH  64.589675   F229B  2023-07-10    F2  29    B    4   5      33.2
1311 194  BH  64.589675   F240B  2023-07-10    F2  40    B    5  12      71.6
1312 194  BH  64.589675   F235D  2023-07-10    F2  35    D    2   1      68.1
1313 194  BH  64.589675    H18B  2023-07-10    H1   8    B    3   9      82.2
1314 194  BH  64.589675   H224B  2023-07-10    H2  24    B    5  29      68.1
1315 194  BH  64.589675   H233B  2023-07-10    H2  33    B    4   7      23.5
1316 194  BH  64.589675   J223D  2023-07-10    J2  23    D    4   9      49.6
1317 194  BH  64.589675   J225D  2023-07-10    J2  25    D    2   5      77.0
1318 194  BH  64.589675   J230C  2023-07-10    J2  30    C    4  14      35.1
1319 194  BH  64.589675   J232D  2023-07-10    J2  32    D    3  12      30.5
1320 194  BH  64.589675   L110A  2023-07-10    L1  10    A    2  14      89.4
1321 194  BH  64.589675    L14D  2023-07-10    L1   4    D    4  10      95.6
1322 194  BH  64.589675   L113C  2023-07-10    L1  13    C    3  13     130.0
1323 194  BH  64.589675   L113D  2023-07-10    L1  13    D    5  22      96.1
1324 194  BH  64.589675   L240D  2023-07-10    L2  40    D    7  13      20.1
1325 194  BH  64.589675   L241D  2023-07-10    L2  41    D    4   1      68.4
1326 194 DPR 112.870361    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1327 195  CC  75.900135    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1328 195  BH  65.126957    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1329 195 DPR 113.062818    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1330 196  CC  75.796951    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1331 196  BH  64.936582    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1332 196 DPR 113.094964    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1333 197  CC  75.792054    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1334 197  BH  64.861278    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1335 197 DPR 113.367115    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1336 198  CC  76.282401    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1337 198  BH  65.414600    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1338 198 DPR 113.767564    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1339 199  CC  76.029065    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1340 199  BH  65.517657    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1341 199 DPR 113.718724    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1342 200  CC  76.450579    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1343 200  BH  65.356284    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1344 200 DPR 113.724729    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1345 201  CC  76.191903    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1346 201  BH  65.556793   D226B  2023-07-17    D2  26    B    5  15      68.2
1347 201  BH  65.556793   D230B  2023-07-17    D2  30    B    4   3      77.8
1348 201  BH  65.556793   D229C  2023-07-17    D2  29    C    4  13      25.7
1349 201  BH  65.556793   D229D  2023-07-17    D2  29    D    2   9      44.2
1350 201  BH  65.556793   D235D  2023-07-17    D2  35    D    3   6      83.6
1351 201  BH  65.556793   F115A  2023-07-17    F1  15    A    6   9     113.2
1352 201  BH  65.556793   F229B  2023-07-17    F2  29    B    4   5      34.3
1353 201  BH  65.556793   F240B  2023-07-17    F2  40    B    5  12      76.3
1354 201  BH  65.556793   F235D  2023-07-17    F2  35    D    2   1      68.8
1355 201  BH  65.556793    H18B  2023-07-17    H1   8    B    3   9      83.7
1356 201  BH  65.556793   H224B  2023-07-17    H2  24    B    5  29      70.3
1357 201  BH  65.556793   H233B  2023-07-17    H2  33    B    4   7      22.9
1358 201  BH  65.556793   J223D  2023-07-17    J2  23    D    4   9      50.8
1359 201  BH  65.556793   J225D  2023-07-17    J2  25    D    2   5      83.5
1360 201  BH  65.556793   J230C  2023-07-17    J2  30    C    4  14      41.6
1361 201  BH  65.556793   J232D  2023-07-17    J2  32    D    3  12      30.7
1362 201  BH  65.556793   L110A  2023-07-17    L1  10    A    2  14      90.0
1363 201  BH  65.556793    L14D  2023-07-17    L1   4    D    4  10      73.2
1364 201  BH  65.556793   L113C  2023-07-17    L1  13    C    3  13     107.6
1365 201  BH  65.556793   L241D  2023-07-17    L2  41    D    4   1      37.4
1366 201 DPR 113.734930    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1367 202  CC  76.114244    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1368 202  BH  65.319181    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1369 202 DPR 114.047121    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1370 203  CC  76.552088    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1371 203  BH  65.345077    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1372 203 DPR 113.826851    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1373 204  CC  76.495529    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1374 204  BH  65.487241    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1375 204 DPR 114.305494    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1376 205  CC  76.298265    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1377 205  BH  65.516874    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1378 205 DPR 113.935431    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1379 206  CC  76.456477    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1380 206  BH  65.474227    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1381 206 DPR 114.302416    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1382 207  CC  76.756729    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1383 207  BH  65.343369    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1384 207 DPR 114.233426    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1385 208  CC  76.390575    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1386 208  BH  65.351869   D226B  2023-07-24    D2  26    B    5  15      68.8
1387 208  BH  65.351869   D230B  2023-07-24    D2  30    B    4   3      79.3
1388 208  BH  65.351869   D229C  2023-07-24    D2  29    C    4  13      26.4
1389 208  BH  65.351869   D229D  2023-07-24    D2  29    D    2   9      44.8
1390 208  BH  65.351869   D235D  2023-07-24    D2  35    D    3   6      86.0
1391 208  BH  65.351869   F115A  2023-07-24    F1  15    A    6   9     109.9
1392 208  BH  65.351869   F229B  2023-07-24    F2  29    B    4   5      36.2
1393 208  BH  65.351869   F240B  2023-07-24    F2  40    B    5  12      75.1
1394 208  BH  65.351869   F235D  2023-07-24    F2  35    D    2   1      69.3
1395 208  BH  65.351869    H18B  2023-07-24    H1   8    B    3   9      81.7
1396 208  BH  65.351869   H224B  2023-07-24    H2  24    B    5  29      68.5
1397 208  BH  65.351869   H233B  2023-07-24    H2  33    B    4   7      22.8
1398 208  BH  65.351869   J225D  2023-07-24    J2  25    D    2   5      81.5
1399 208  BH  65.351869   J230C  2023-07-24    J2  30    C    4  14      32.8
1400 208  BH  65.351869   J232D  2023-07-24    J2  32    D    3  12      30.2
1401 208  BH  65.351869   L110A  2023-07-24    L1  10    A    2  14      91.2
1402 208  BH  65.351869    L14D  2023-07-24    L1   4    D    4  10      93.3
1403 208  BH  65.351869   L113C  2023-07-24    L1  13    C    3  13      89.8
1404 208  BH  65.351869   L113D  2023-07-24    L1  13    D    5  22     112.1
1405 208  BH  65.351869   L241D  2023-07-24    L2  41    D    4   1      69.6
1406 208 DPR 114.375749    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
     herbivory     pred2
1         <NA>  3.967930
2         <NA>  3.992824
3         <NA>  4.237680
4            Y  4.201579
5         <NA>  4.430940
6            Y  3.824152
7            Y  3.982101
8         <NA>  3.975611
9         <NA>  3.841476
10        <NA>  3.805019
11           Y  3.969249
12        <NA>  3.972846
13           Y  3.710265
14        <NA>  3.803876
15        <NA>  4.209246
16        <NA>  4.471308
17           Y  4.186500
18        <NA>  3.778138
19           Y  3.932197
20           Y  4.347765
21        <NA>  4.127186
22        <NA>  4.114924
23        <NA>  4.038419
24        <NA>  3.922018
25        <NA>  3.963204
26        <NA>  3.888205
27        <NA>  3.975314
28           Y  4.233056
29        <NA>  4.060493
30        <NA>  3.742677
31           Y  4.201858
32        <NA>  4.303181
33        <NA>  3.852062
34        <NA>  4.339974
35        <NA>  4.324464
36           Y  4.132320
37        <NA>  4.165742
38        <NA>  3.960452
39        <NA>  4.096894
40        <NA>  4.177245
41        <NA>  4.214643
42        <NA>  4.099426
43        <NA>  3.936184
44        <NA>  3.811925
45        <NA>  4.117418
46        <NA>  3.806436
47        <NA>  4.033292
48        <NA>  3.914814
49        <NA>  4.167884
50        <NA>  4.173428
51           Y  4.183971
52        <NA>  4.285362
53        <NA>  4.090425
54        <NA>  4.346586
55        <NA>  4.194884
56           Y  3.906552
57        <NA>  4.150396
58        <NA>  3.656004
59        <NA>  4.025263
60        <NA>  3.869618
61           Y  4.162538
62        <NA>  4.282597
63        <NA>  3.946302
64        <NA>  4.112193
65        <NA>  3.685549
66        <NA>  3.686741
67        <NA>  4.047531
68        <NA>  3.930772
69        <NA>  3.923525
70           Y  4.069708
71        <NA>  4.117482
72        <NA>  3.877490
73           Y  4.146556
74        <NA>  4.296136
75           Y  4.016954
76           Y  4.340213
77           Y  3.970090
78           Y  4.261809
79        <NA>  3.992450
80           Y  3.758761
81        <NA>  3.984880
82        <NA>  3.859404
83        <NA>  3.975836
84           Y  3.862840
85           Y  4.253840
86           Y  3.860294
87        <NA>  3.718866
88        <NA>  3.749390
89           Y  3.739358
90        <NA>  4.104320
91           Y  4.122123
92           Y  4.009375
93           Y  4.240476
94           Y  3.740548
95        <NA>  3.833363
96           Y  4.176585
97           Y  4.246419
98        <NA>  3.842487
99        <NA>  3.891785
100          Y  3.867265
101       <NA>  3.990957
102          Y  4.418764
103       <NA>  4.088692
104       <NA>  4.047897
105       <NA>  4.253043
106          Y  3.974667
107       <NA>  3.915674
108       <NA>  3.793520
109          Y  4.166034
110       <NA>  3.745286
111       <NA>  4.194429
112       <NA>  4.162046
113       <NA>  4.117195
114          Y  4.220344
115          Y  3.751292
116       <NA>  3.928635
117       <NA>  3.820682
118       <NA>  4.098745
119          Y  3.923083
120       <NA>  4.136613
121          Y  4.060423
122       <NA>  3.864438
123       <NA>  3.755844
124       <NA>  3.918103
125       <NA>  3.907934
126          Y  3.798792
127          Y  4.064719
128          Y  3.996677
129       <NA>  3.976673
130          Y  3.897686
131       <NA>  3.922574
132          Y  4.152092
133          Y  4.031853
134       <NA>  4.113065
135          Y  4.202572
136       <NA>  4.190770
137          Y  4.044949
138          Y  4.078987
139       <NA>  4.079329
140       <NA>  3.843873
141          Y  3.896594
142       <NA>  3.952706
143       <NA>  4.004008
144       <NA>  3.945877
145          Y  4.006237
146       <NA>  3.795776
147          Y  4.061952
148       <NA>  4.026067
149       <NA>  3.978190
150       <NA>  4.078230
151       <NA>  3.843217
152       <NA>  4.119677
153       <NA>  4.008003
154       <NA>  3.991561
155          Y  4.040298
156       <NA>        NA
157       <NA>        NA
158       <NA>        NA
159       <NA>        NA
160       <NA>        NA
161       <NA>        NA
162       <NA>        NA
163       <NA>        NA
164       <NA>        NA
165       <NA>        NA
166       <NA>        NA
167       <NA>        NA
168       <NA>        NA
169       <NA>        NA
170       <NA>        NA
171       <NA>        NA
172       <NA>        NA
173       <NA>        NA
174       <NA>        NA
175       <NA>        NA
176       <NA>        NA
177       <NA>        NA
178       <NA>        NA
179       <NA>        NA
180       <NA>        NA
181       <NA>        NA
182       <NA>        NA
183       <NA>        NA
184       <NA>        NA
185       <NA>        NA
186       <NA>        NA
187       <NA>        NA
188       <NA>        NA
189       <NA>        NA
190       <NA>        NA
191       <NA>        NA
192       <NA>        NA
193       <NA>        NA
194       <NA>        NA
195       <NA>        NA
196       <NA>        NA
197       <NA>        NA
198       <NA>        NA
199       <NA>        NA
200       <NA>        NA
201       <NA>        NA
202       <NA>        NA
203       <NA>        NA
204       <NA>        NA
205       <NA>        NA
206       <NA>        NA
207       <NA>        NA
208       <NA>        NA
209       <NA>        NA
210       <NA>        NA
211       <NA>        NA
212       <NA>        NA
213       <NA>        NA
214       <NA>        NA
215       <NA>        NA
216       <NA>  3.956034
217       <NA>  4.250841
218       <NA>  4.206485
219       <NA>  3.849379
220       <NA>  4.093494
221       <NA>  4.222301
222       <NA>  4.226680
223       <NA>  3.818979
224       <NA>  4.039906
225       <NA>  4.139798
226       <NA>  4.295933
227       <NA>  3.812357
228       <NA>  4.286890
229       <NA>  4.096165
230       <NA>  4.584393
231       <NA>  4.288621
232       <NA>  4.355224
233          N  3.711352
234       <NA>  3.965897
235       <NA>  4.250089
236       <NA>  4.069742
237       <NA>  4.246811
238       <NA>  4.153869
239       <NA>  4.048554
240       <NA>  4.136725
241       <NA>  4.273495
242       <NA>  4.235323
243       <NA>  4.339068
244       <NA>  4.141982
245       <NA>  3.994703
246       <NA>  4.149279
247       <NA>  3.961219
248       <NA>  4.377505
249       <NA>  4.228358
250       <NA>  4.242585
251       <NA>  4.038363
252       <NA>  4.067020
253       <NA>  3.815169
254       <NA>  3.910846
255       <NA>  4.203260
256       <NA>  3.958338
257       <NA>  4.439005
258       <NA>  4.264494
259       <NA>  3.855434
260       <NA>  3.855946
261       <NA>  3.992318
262       <NA>  4.356212
263       <NA>  3.951115
264       <NA>  3.928985
265       <NA>  3.976030
266       <NA>  3.969985
267       <NA>  4.091435
268       <NA>  4.081862
269       <NA>  4.122894
270       <NA>  4.587742
271       <NA>  3.879176
272       <NA>  3.762552
273       <NA>  4.014288
274          N  3.922272
275       <NA>  4.127760
276       <NA>  3.946959
277       <NA>  4.211549
278       <NA>  4.113569
279       <NA>  3.975225
280       <NA>  3.924323
281       <NA>  4.144186
282       <NA>  3.886072
283       <NA>  3.851816
284       <NA>  4.023414
285       <NA>  4.151857
286       <NA>  4.227586
287       <NA>  3.783616
288       <NA>  4.168561
289       <NA>  4.056065
290       <NA>  4.060807
291       <NA>  4.393031
292       <NA>  4.112142
293       <NA>  4.389499
294       <NA>  4.243682
295       <NA>  4.168721
296       <NA>  4.178240
297       <NA>  4.379732
298       <NA>  4.139031
299       <NA>  4.171182
300       <NA>  4.152611
301       <NA>  4.092606
302       <NA>  4.180620
303       <NA>  3.921985
304       <NA>  4.132516
305       <NA>  4.043901
306       <NA>  4.126033
307       <NA>  4.152787
308       <NA>  4.153038
309       <NA>  4.149765
310       <NA>  4.547556
311       <NA>  3.968068
312       <NA>  4.233018
313       <NA>  4.036543
314       <NA>  4.155202
315       <NA>  3.962074
316       <NA>  4.210530
317       <NA>  4.318203
318       <NA>  4.395005
319       <NA>  4.127880
320       <NA>  4.347502
321       <NA>  3.742864
322       <NA>  3.935772
323       <NA>  4.114292
324       <NA>  4.050114
325       <NA>  3.943647
326       <NA>  4.179904
327       <NA>  3.911197
328       <NA>  4.245165
329       <NA>  3.937842
330       <NA>  4.069867
331       <NA>  4.082403
332       <NA>  4.542889
333       <NA>  4.383129
334       <NA>  3.962815
335       <NA>  4.152239
336       <NA>  4.174636
337       <NA>  4.401785
338       <NA>  4.294124
339       <NA>  4.442232
340       <NA>  4.080530
341       <NA>  4.372652
342          2  4.014971
343       <NA>        NA
344       <NA>        NA
345       <NA>        NA
346       <NA>        NA
347       <NA>        NA
348       <NA>        NA
349       <NA>        NA
350       <NA>        NA
351       <NA>        NA
352       <NA>        NA
353       <NA>        NA
354       <NA>        NA
355       <NA>        NA
356       <NA>        NA
357       <NA>        NA
358       <NA>        NA
359       <NA>        NA
360       <NA>        NA
361       <NA>        NA
362       <NA>        NA
363       <NA>        NA
364       <NA>        NA
365       <NA>        NA
366       <NA>        NA
367       <NA>        NA
368       <NA>        NA
369       <NA>        NA
370       <NA>        NA
371       <NA>        NA
372       <NA>        NA
373       <NA>        NA
374       <NA>        NA
375       <NA>        NA
376       <NA>        NA
377       <NA>        NA
378       <NA>        NA
379       <NA>        NA
380       <NA>        NA
381       <NA>        NA
382       <NA>  4.645039
383       <NA>  4.650564
384       <NA>  4.788565
385       <NA>  4.071687
386       <NA>  4.754207
387       <NA>  4.421756
388       <NA>  4.255781
389       <NA>  4.548027
390       <NA>  4.639866
391       <NA>  4.534559
392       <NA>  4.407716
393       <NA>  4.592904
394       <NA>  4.462216
395       <NA>  4.788584
396       <NA>  4.378211
397       <NA>  4.751840
398       <NA>  4.271627
399       <NA>  4.474034
400       <NA>  4.340750
401       <NA>  4.324762
402       <NA>  4.467728
403       <NA>  4.659763
404       <NA>  4.553852
405       <NA>  4.383504
406       <NA>  4.563212
407       <NA>  4.437894
408       <NA>  4.850357
409       <NA>  4.191552
410       <NA>  4.385826
411       <NA>  4.447791
412       <NA>  4.317379
413       <NA>  4.459160
414       <NA>  4.224985
415       <NA>  4.474080
416       <NA>  4.219580
417       <NA>  4.467479
418       <NA>  4.501575
419       <NA>  4.519510
420       <NA>  4.491163
421       <NA>  4.489510
422       <NA>  4.552617
423       <NA>  4.594708
424       <NA>  4.238385
425       <NA>  4.434393
426       <NA>  4.631661
427       <NA>  4.072933
428       <NA>  4.361179
429       <NA>  4.468451
430       <NA>  4.259791
431       <NA>  4.483213
432       <NA>  4.170975
433       <NA>  3.926828
434       <NA>  4.328804
435       <NA>  4.302769
436       <NA>  4.188455
437       <NA>  4.216736
438       <NA>  4.601528
439       <NA>  3.984222
440       <NA>  4.476863
441       <NA>  4.429721
442       <NA>  4.511176
443       <NA>  4.469856
444       <NA>  4.540613
445       <NA>  4.471158
446       <NA>  4.387375
447       <NA>  4.062076
448       <NA>  4.236269
449       <NA>  4.573770
450       <NA>  4.254848
451       <NA>  4.561226
452       <NA>  4.656360
453       <NA>  4.586994
454       <NA>  4.438618
455       <NA>  4.111678
456       <NA>  4.615214
457       <NA>  4.378463
458       <NA>  4.070121
459       <NA>  4.436832
460       <NA>  4.406292
461       <NA>  4.526834
462       <NA>  4.704886
463       <NA>  4.579209
464       <NA>  4.456807
465       <NA>  4.818778
466       <NA>  4.862850
467       <NA>  4.624944
468       <NA>  4.506260
469       <NA>  5.017193
470       <NA>  4.882535
471       <NA>  4.913032
472       <NA>  4.647868
473       <NA>  4.545386
474       <NA>  4.858841
475       <NA>  4.892642
476       <NA>        NA
477       <NA>        NA
478       <NA>        NA
479       <NA>        NA
480       <NA>        NA
481       <NA>        NA
482       <NA>        NA
483       <NA>        NA
484       <NA>        NA
485       <NA>        NA
486       <NA>        NA
487       <NA>        NA
488       <NA>        NA
489       <NA>        NA
490       <NA>        NA
491       <NA>        NA
492       <NA>        NA
493       <NA>        NA
494       <NA>        NA
495       <NA>        NA
496       <NA>        NA
497       <NA>        NA
498       <NA>        NA
499       <NA>        NA
500       <NA>        NA
501       <NA>        NA
502       <NA>        NA
503       <NA>        NA
504       <NA>        NA
505       <NA>        NA
506       <NA>        NA
507       <NA>        NA
508       <NA>        NA
509       <NA>        NA
510       <NA>        NA
511       <NA>        NA
512       <NA>        NA
513       <NA>        NA
514       <NA>        NA
515       <NA>  5.157322
516       <NA>  5.601636
517       <NA>  5.610458
518       <NA>  5.450774
519       <NA>  5.579050
520       <NA>  5.497488
521       <NA>  5.227321
522       <NA>  5.477058
523       <NA>  5.177994
524       <NA>  5.567473
525       <NA>  5.661558
526       <NA>  5.589363
527       <NA>  5.570243
528       <NA>  5.573044
529       <NA>  5.086037
530       <NA>  4.980489
531       <NA>  5.604455
532       <NA>  5.114610
533       <NA>  5.172907
534       <NA>  5.289276
535       <NA>  5.422354
536       <NA>  5.320595
537       <NA>  5.088948
538       <NA>  5.462119
539       <NA>  5.313172
540       <NA>  5.572975
541       <NA>  5.346363
542       <NA>  5.571092
543       <NA>  5.527262
544       <NA>  5.540419
545       <NA>  5.356987
546       <NA>  5.125891
547       <NA>  5.409279
548       <NA>  5.146534
549       <NA>  5.323113
550       <NA>  5.588955
551       <NA>  5.218473
552       <NA>  5.346870
553       <NA>  5.085513
554       <NA>  5.479904
555       <NA>  5.136637
556       <NA>  5.267807
557       <NA>  5.232099
558       <NA>  5.441368
559       <NA>  5.389647
560       <NA>  5.298034
561       <NA>  5.401664
562       <NA>  5.188026
563       <NA>  5.285414
564       <NA>  5.272613
565       <NA>  5.035810
566       <NA>  5.293219
567       <NA>  5.080928
568       <NA>  5.346249
569       <NA>  5.487651
570       <NA>  5.068518
571       <NA>  5.478759
572       <NA>  5.029995
573       <NA>  5.557400
574       <NA>  5.302079
575       <NA>  4.947259
576       <NA>  5.402394
577       <NA>  5.249170
578       <NA>  5.606184
579       <NA>  5.299818
580       <NA>  4.967197
581       <NA>  5.267040
582       <NA>  5.596362
583       <NA>  4.922482
584       <NA>  6.504796
585       <NA>  6.128811
586       <NA>  6.419585
587       <NA>  6.225083
588       <NA>  6.572256
589       <NA>  6.322072
590       <NA>  6.176444
591       <NA>  6.214700
592       <NA>  6.276687
593       <NA>        NA
594       <NA>        NA
595       <NA>        NA
596       <NA>        NA
597       <NA>        NA
598       <NA>        NA
599       <NA>        NA
600       <NA>        NA
601       <NA>        NA
602       <NA>        NA
603       <NA>        NA
604       <NA>        NA
605       <NA>        NA
606       <NA>        NA
607       <NA>        NA
608       <NA>        NA
609       <NA>        NA
610       <NA>        NA
611          N  6.485304
612       <NA>  6.307636
613       <NA>  6.408489
614       <NA>  6.790899
615       <NA>  6.414138
616          Y  6.438446
617       <NA>  6.456467
618       <NA>  6.478001
619          N  6.396759
620          N  6.428626
621       <NA>  6.527634
622       <NA>  5.714984
623       <NA>  6.100314
624       <NA>  5.931933
625       <NA>  6.038090
626       <NA>  6.104848
627       <NA>  6.082096
628       <NA>  6.236007
629          Y  5.925827
630       <NA>  6.231060
631       <NA>  5.991361
632       <NA>  6.164238
633       <NA>  5.828256
634       <NA>  6.266452
635       <NA>  5.977702
636       <NA>  6.091193
637       <NA>  6.024175
638       <NA>  5.917834
639       <NA>  6.061683
640          Y  6.258178
641          Y  6.206826
642       <NA>  6.159607
643          Y  5.939701
644       <NA>  5.666478
645       <NA>  6.272509
646       <NA>  5.956701
647       <NA>  6.050409
648       <NA>  6.223828
649       <NA>  6.044855
650       <NA>  5.942472
651       <NA>  5.961268
652       <NA>  6.232720
653       <NA>  5.866648
654       <NA>  6.091708
655       <NA>  5.843695
656       <NA>  6.218626
657       <NA>  5.892235
658       <NA>  6.061425
659       <NA>  5.768552
660       <NA>  6.009130
661       <NA>  6.060646
662       <NA>  5.851265
663       <NA>  6.618512
664          Y  5.910208
665       <NA>  5.854919
666       <NA>  5.984828
667       <NA>  7.634706
668       <NA>  7.524079
669       <NA>  7.460225
670       <NA>  7.733761
671       <NA>  7.777300
672       <NA>  7.673781
673       <NA>  7.925711
674          N  7.777242
675       <NA>        NA
676       <NA>        NA
677       <NA>        NA
678       <NA>        NA
679       <NA>        NA
680       <NA>        NA
681       <NA>        NA
682       <NA>        NA
683       <NA>        NA
684       <NA>        NA
685       <NA>        NA
686       <NA>        NA
687       <NA>        NA
688       <NA>        NA
689       <NA>        NA
690       <NA>        NA
691       <NA>        NA
692       <NA>        NA
693       <NA>        NA
694       <NA>        NA
695       <NA>        NA
696       <NA>        NA
697       <NA>        NA
698       <NA>        NA
699       <NA>        NA
700       <NA>        NA
701       <NA>        NA
702       <NA>  8.741736
703       <NA>  8.353367
704       <NA>  8.723756
705       <NA>  8.336899
706       <NA>  8.729912
707       <NA>  8.546976
708       <NA>  7.808558
709       <NA>  7.822508
710       <NA>  7.860971
711       <NA>  8.002398
712       <NA>  7.525036
713       <NA>  7.773971
714       <NA>  8.423394
715       <NA>  7.935282
716       <NA>  8.067286
717       <NA>  7.720240
718       <NA>  7.942433
719       <NA>  7.908721
720       <NA>  7.939602
721       <NA>  7.874263
722       <NA>  7.898484
723       <NA>  7.924880
724       <NA>  8.034522
725       <NA>  7.547507
726       <NA>  7.878305
727       <NA>  7.809313
728       <NA>  7.756488
729       <NA>  7.997182
730       <NA>  8.268266
731       <NA>  7.919449
732       <NA>  7.708223
733       <NA>  8.066466
734       <NA>  8.077956
735       <NA>  7.800708
736       <NA>  7.811560
737       <NA>  7.689406
738       <NA>  8.129433
739       <NA>  8.026164
740       <NA>  7.914904
741       <NA>  7.781555
742       <NA>  8.029073
743       <NA>  8.180117
744       <NA>  8.158895
745       <NA>  8.047891
746       <NA>  7.537409
747       <NA>  8.060977
748       <NA>  7.808091
749       <NA>  7.822330
750       <NA> 10.987032
751       <NA> 11.127627
752       <NA> 11.350289
753          Y 11.323896
754       <NA> 10.914931
755       <NA>        NA
756       <NA>        NA
757       <NA>        NA
758       <NA>        NA
759       <NA>        NA
760       <NA>        NA
761       <NA>        NA
762       <NA>        NA
763       <NA>        NA
764       <NA>        NA
765       <NA>        NA
766       <NA>        NA
767       <NA>        NA
768       <NA>        NA
769       <NA>        NA
770       <NA>        NA
771       <NA>        NA
772       <NA>        NA
773       <NA> 11.045267
774       <NA> 10.930731
775       <NA> 10.756795
776       <NA> 10.538185
777       <NA> 10.962599
778       <NA>  9.941635
779       <NA> 10.031908
780       <NA>  9.783300
781       <NA>  9.744320
782          Y  9.873110
783       <NA>  9.780779
784       <NA>  9.770609
785       <NA> 10.039914
786       <NA> 10.263066
787       <NA>  9.911083
788       <NA>  9.821457
789       <NA>  9.949486
790       <NA>  9.649868
791       <NA> 10.052981
792       <NA>  9.811011
793       <NA>  9.882549
794       <NA>  9.870599
795          Y  9.967807
796       <NA>  9.609889
797       <NA>  9.934134
798       <NA> 10.066971
799       <NA>  9.820105
800       <NA>  9.523714
801       <NA>  9.899188
802       <NA>  9.639840
803       <NA> 10.110614
804       <NA>  9.959017
805       <NA>  9.747875
806       <NA>  9.998054
807       <NA>  9.857868
808       <NA> 10.213880
809       <NA>  9.779760
810       <NA>  9.912668
811       <NA>  9.779674
812       <NA>  9.864551
813       <NA> 10.043419
814       <NA> 10.000832
815       <NA>  9.690564
816       <NA>  9.669143
817       <NA> 14.603094
818       <NA> 14.634118
819       <NA> 14.453849
820       <NA> 14.679684
821       <NA>        NA
822       <NA>        NA
823       <NA>        NA
824       <NA>        NA
825       <NA>        NA
826       <NA>        NA
827       <NA>        NA
828       <NA>        NA
829       <NA>        NA
830       <NA>        NA
831       <NA>        NA
832       <NA>        NA
833       <NA>        NA
834       <NA>        NA
835       <NA>        NA
836       <NA>        NA
837       <NA>        NA
838       <NA>        NA
839       <NA> 14.115381
840       <NA> 13.918731
841       <NA> 14.208390
842       <NA> 13.962165
843       <NA> 12.581319
844       <NA> 12.540335
845       <NA> 12.553030
846       <NA> 12.351847
847       <NA> 12.002117
848       <NA> 12.380310
849       <NA> 12.374347
850       <NA> 12.632744
851       <NA> 12.382552
852       <NA> 12.309866
853       <NA> 12.726592
854       <NA> 12.729922
855       <NA> 12.135058
856       <NA> 12.325014
857       <NA> 12.578090
858       <NA> 12.526711
859       <NA> 12.711437
860       <NA> 12.589215
861       <NA> 12.555437
862       <NA> 12.409455
863       <NA> 12.222346
864       <NA> 12.572187
865       <NA> 12.542645
866       <NA> 12.516587
867       <NA> 12.453451
868       <NA> 12.532166
869       <NA> 12.479528
870       <NA> 12.426055
871       <NA> 12.449242
872       <NA> 12.702289
873       <NA> 12.469210
874       <NA> 12.480333
875       <NA> 12.455870
876       <NA> 12.483776
877       <NA> 12.724485
878       <NA> 12.719261
879       <NA> 12.273399
880       <NA> 19.430770
881       <NA> 19.335687
882       <NA> 19.337160
883       <NA> 19.256395
884       <NA>        NA
885       <NA>        NA
886       <NA>        NA
887       <NA>        NA
888       <NA>        NA
889       <NA>        NA
890       <NA>        NA
891       <NA>        NA
892       <NA>        NA
893       <NA>        NA
894       <NA>        NA
895       <NA>        NA
896       <NA>        NA
897       <NA>        NA
898       <NA>        NA
899       <NA>        NA
900       <NA>        NA
901       <NA>        NA
902       <NA> 18.253644
903       <NA> 17.804625
904       <NA> 17.767766
905       <NA> 15.662985
906       <NA> 15.844135
907          N 15.749529
908       <NA> 15.869989
909       <NA> 15.900161
910       <NA> 15.981030
911       <NA> 15.945625
912       <NA> 16.048468
913       <NA> 15.941699
914       <NA> 15.694985
915       <NA> 15.894896
916       <NA> 15.938810
917       <NA> 16.034355
918       <NA> 15.902064
919       <NA> 15.853825
920       <NA> 16.000820
921       <NA> 16.105806
922       <NA> 15.927681
923       <NA> 15.963269
924       <NA> 15.591900
925       <NA> 16.084805
926       <NA> 15.874943
927       <NA> 16.184136
928       <NA> 15.729415
929       <NA> 15.925902
930       <NA> 16.095084
931       <NA> 15.985796
932       <NA> 16.108695
933       <NA> 16.064907
934       <NA> 16.294609
935       <NA> 15.884684
936       <NA> 15.767750
937          N 15.688354
938       <NA> 16.093076
939       <NA> 15.882409
940       <NA> 16.163009
941       <NA> 25.652379
942       <NA> 25.456939
943       <NA> 25.166923
944       <NA> 25.632342
945       <NA>        NA
946       <NA>        NA
947       <NA>        NA
948       <NA>        NA
949       <NA>        NA
950       <NA>        NA
951       <NA>        NA
952       <NA>        NA
953       <NA>        NA
954       <NA>        NA
955       <NA>        NA
956       <NA>        NA
957       <NA>        NA
958       <NA>        NA
959       <NA>        NA
960       <NA>        NA
961       <NA>        NA
962       <NA>        NA
963       <NA> 22.850965
964       <NA> 23.028561
965       <NA> 23.171361
966       <NA> 20.202022
967       <NA> 20.318193
968       <NA> 20.192909
969       <NA> 20.104512
970       <NA> 20.226745
971          Y 20.322474
972       <NA> 20.394543
973       <NA> 20.135689
974       <NA> 20.408189
975       <NA> 20.410134
976       <NA> 20.178404
977       <NA> 20.438999
978       <NA> 20.279600
979          Y 19.800695
980       <NA> 20.236882
981       <NA> 20.515095
982       <NA> 20.153181
983       <NA> 20.222603
984       <NA> 20.051768
985       <NA> 20.119145
986       <NA> 20.126332
987       <NA> 20.248344
988       <NA> 20.332439
989       <NA> 20.123860
990       <NA> 20.072000
991       <NA> 20.260424
992       <NA> 20.311991
993       <NA> 20.068870
994       <NA> 20.460420
995       <NA> 20.214084
996       <NA> 20.239446
997       <NA> 20.310543
998       <NA> 20.188166
999       <NA> 20.468886
1000      <NA> 19.926164
1001      <NA> 20.096924
1002      <NA> 33.166234
1003      <NA> 32.911488
1004      <NA> 33.140779
1005      <NA>        NA
1006      <NA>        NA
1007      <NA>        NA
1008      <NA>        NA
1009      <NA>        NA
1010      <NA>        NA
1011      <NA>        NA
1012      <NA>        NA
1013      <NA>        NA
1014      <NA>        NA
1015      <NA>        NA
1016      <NA>        NA
1017      <NA>        NA
1018      <NA>        NA
1019      <NA>        NA
1020      <NA>        NA
1021      <NA>        NA
1022      <NA>        NA
1023      <NA> 29.758816
1024      <NA> 25.452562
1025      <NA> 25.290039
1026      <NA> 25.431052
1027      <NA> 25.717934
1028      <NA> 25.816319
1029      <NA> 25.449219
1030      <NA> 25.339557
1031      <NA> 25.065208
1032      <NA> 25.297956
1033      <NA> 25.486261
1034      <NA> 25.409989
1035      <NA> 25.610614
1036      <NA> 25.159067
1037      <NA> 25.490181
1038      <NA> 25.573482
1039      <NA> 25.179008
1040      <NA> 25.411708
1041      <NA> 25.415644
1042      <NA> 25.117641
1043      <NA> 25.625072
1044      <NA> 25.485259
1045      <NA> 25.386334
1046      <NA> 25.087733
1047      <NA> 25.342553
1048      <NA> 25.722929
1049      <NA> 25.210836
1050      <NA> 25.507118
1051      <NA> 25.625197
1052      <NA> 25.287765
1053      <NA> 25.428525
1054      <NA> 25.152489
1055      <NA> 25.514673
1056      <NA> 25.190733
1057      <NA> 25.421622
1058      <NA> 25.356321
1059      <NA> 25.553842
1060      <NA> 42.501069
1061      <NA>        NA
1062      <NA>        NA
1063      <NA>        NA
1064      <NA>        NA
1065      <NA>        NA
1066      <NA>        NA
1067      <NA>        NA
1068      <NA>        NA
1069      <NA>        NA
1070      <NA>        NA
1071      <NA>        NA
1072      <NA>        NA
1073      <NA>        NA
1074      <NA>        NA
1075      <NA>        NA
1076      <NA>        NA
1077      <NA>        NA
1078      <NA>        NA
1079      <NA>        NA
1080      <NA>        NA
1081      <NA>        NA
1082      <NA>        NA
1083      <NA>        NA
1084      <NA>        NA
1085      <NA>        NA
1086      <NA>        NA
1087      <NA>        NA
1088      <NA>        NA
1089      <NA>        NA
1090      <NA>        NA
1091      <NA>        NA
1092      <NA>        NA
1093      <NA>        NA
1094      <NA>        NA
1095      <NA>        NA
1096      <NA>        NA
1097      <NA>        NA
1098      <NA>        NA
1099      <NA>        NA
1100      <NA>        NA
1101      <NA>        NA
1102      <NA>        NA
1103      <NA>        NA
1104      <NA>        NA
1105      <NA>        NA
1106      <NA>        NA
1107      <NA>        NA
1108      <NA>        NA
1109      <NA>        NA
1110      <NA>        NA
1111      <NA>        NA
1112      <NA>        NA
1113      <NA>        NA
1114      <NA>        NA
1115      <NA>        NA
1116      <NA>        NA
1117      <NA>        NA
1118      <NA>        NA
1119      <NA>        NA
1120      <NA>        NA
1121      <NA>        NA
1122      <NA>        NA
1123      <NA>        NA
1124      <NA>        NA
1125      <NA>        NA
1126      <NA>        NA
1127      <NA>        NA
1128      <NA>        NA
1129      <NA>        NA
1130      <NA>        NA
1131      <NA>        NA
1132      <NA>        NA
1133      <NA>        NA
1134      <NA>        NA
1135      <NA>        NA
1136      <NA>        NA
1137      <NA>        NA
1138      <NA>        NA
1139      <NA>        NA
1140      <NA>        NA
1141      <NA>        NA
1142      <NA>        NA
1143      <NA> 49.611186
1144      <NA> 50.090009
1145      <NA> 49.668964
1146      <NA> 50.034470
1147      <NA> 49.772706
1148      <NA> 50.226468
1149      <NA> 49.637030
1150      <NA> 49.734929
1151      <NA> 49.641693
1152      <NA> 49.573821
1153      <NA> 49.823789
1154      <NA> 49.462741
1155      <NA> 49.472238
1156      <NA> 49.660685
1157      <NA> 49.786765
1158      <NA> 49.709176
1159      <NA> 49.570911
1160      <NA> 49.591908
1161      <NA> 50.117569
1162      <NA> 49.871700
1163      <NA> 49.693434
1164      <NA> 49.497883
1165      <NA> 49.988336
1166      <NA> 49.824434
1167      <NA> 50.255251
1168      <NA> 49.745790
1169      <NA> 49.767684
1170      <NA> 50.071958
1171      <NA> 49.783756
1172      <NA> 49.652292
1173      <NA>        NA
1174      <NA>        NA
1175      <NA>        NA
1176      <NA>        NA
1177      <NA>        NA
1178      <NA>        NA
1179      <NA>        NA
1180      <NA>        NA
1181      <NA>        NA
1182      <NA>        NA
1183      <NA>        NA
1184      <NA>        NA
1185      <NA>        NA
1186      <NA>        NA
1187      <NA>        NA
1188      <NA>        NA
1189      <NA>        NA
1190      <NA>        NA
1191      <NA>        NA
1192      <NA>        NA
1193      <NA>        NA
1194      <NA>        NA
1195      <NA>        NA
1196      <NA>        NA
1197      <NA>        NA
1198      <NA>        NA
1199      <NA>        NA
1200      <NA>        NA
1201      <NA>        NA
1202      <NA>        NA
1203      <NA>        NA
1204      <NA>        NA
1205      <NA>        NA
1206      <NA>        NA
1207      <NA>        NA
1208      <NA>        NA
1209      <NA>        NA
1210      <NA>        NA
1211      <NA>        NA
1212      <NA>        NA
1213      <NA>        NA
1214      <NA>        NA
1215      <NA>        NA
1216      <NA>        NA
1217      <NA>        NA
1218      <NA>        NA
1219      <NA>        NA
1220      <NA>        NA
1221      <NA>        NA
1222      <NA>        NA
1223      <NA>        NA
1224      <NA>        NA
1225      <NA>        NA
1226      <NA>        NA
1227      <NA>        NA
1228      <NA>        NA
1229      <NA>        NA
1230      <NA>        NA
1231      <NA>        NA
1232      <NA>        NA
1233      <NA>        NA
1234      <NA>        NA
1235      <NA>        NA
1236      <NA>        NA
1237      <NA>        NA
1238      <NA>        NA
1239      <NA>        NA
1240      <NA>        NA
1241      <NA>        NA
1242      <NA>        NA
1243      <NA>        NA
1244      <NA>        NA
1245      <NA>        NA
1246      <NA>        NA
1247      <NA>        NA
1248      <NA>        NA
1249      <NA>        NA
1250      <NA>        NA
1251      <NA>        NA
1252      <NA>        NA
1253      <NA>        NA
1254      <NA>        NA
1255      <NA>        NA
1256      <NA> 63.819263
1257      <NA> 63.616593
1258      <NA> 63.997405
1259      <NA> 63.912174
1260      <NA> 63.567954
1261      <NA> 63.690182
1262      <NA> 63.804713
1263      <NA> 63.439691
1264      <NA> 63.513517
1265      <NA> 63.899540
1266      <NA> 63.557876
1267      <NA> 63.617638
1268      <NA> 63.952739
1269      <NA> 64.237896
1270      <NA> 63.690648
1271      <NA> 63.590982
1272      <NA> 63.789974
1273      <NA> 63.667381
1274      <NA> 63.472919
1275      <NA> 63.668144
1276      <NA> 63.511648
1277      <NA> 63.747202
1278      <NA> 63.958576
1279      <NA> 63.839727
1280      <NA> 63.745304
1281      <NA> 63.358017
1282      <NA> 63.581919
1283      <NA> 63.068750
1284      <NA>        NA
1285      <NA>        NA
1286      <NA>        NA
1287      <NA>        NA
1288      <NA>        NA
1289      <NA>        NA
1290      <NA>        NA
1291      <NA>        NA
1292      <NA>        NA
1293      <NA>        NA
1294      <NA>        NA
1295      <NA>        NA
1296      <NA>        NA
1297      <NA>        NA
1298      <NA>        NA
1299      <NA>        NA
1300      <NA>        NA
1301      <NA>        NA
1302      <NA>        NA
1303      <NA>        NA
1304      <NA> 64.749143
1305      <NA> 64.759712
1306      <NA> 64.551668
1307      <NA> 64.669087
1308      <NA> 64.699616
1309      <NA> 64.538542
1310      <NA> 64.531808
1311      <NA> 64.934681
1312      <NA> 64.813888
1313      <NA> 64.629411
1314      <NA> 64.933353
1315      <NA> 64.520234
1316      <NA> 64.964933
1317      <NA> 64.797395
1318      <NA> 64.423291
1319      <NA> 64.343010
1320      <NA> 64.674172
1321      <NA> 64.902631
1322      <NA> 64.669107
1323      <NA> 64.454531
1324      <NA> 64.737150
1325      <NA> 64.436968
1326      <NA>        NA
1327      <NA>        NA
1328      <NA>        NA
1329      <NA>        NA
1330      <NA>        NA
1331      <NA>        NA
1332      <NA>        NA
1333      <NA>        NA
1334      <NA>        NA
1335      <NA>        NA
1336      <NA>        NA
1337      <NA>        NA
1338      <NA>        NA
1339      <NA>        NA
1340      <NA>        NA
1341      <NA>        NA
1342      <NA>        NA
1343      <NA>        NA
1344      <NA>        NA
1345      <NA>        NA
1346      <NA> 65.064142
1347      <NA> 65.046083
1348      <NA> 65.324325
1349      <NA> 65.021228
1350      <NA> 65.652136
1351      <NA> 65.186478
1352      <NA> 64.843652
1353      <NA> 65.391702
1354      <NA> 65.252309
1355      <NA> 65.329295
1356      <NA> 65.112609
1357      <NA> 65.379202
1358      <NA> 65.030391
1359      <NA> 65.216103
1360      <NA> 65.433353
1361      <NA> 65.435194
1362      <NA> 65.308149
1363      <NA> 65.060312
1364      <NA> 65.210076
1365      <NA> 65.103566
1366      <NA>        NA
1367      <NA>        NA
1368      <NA>        NA
1369      <NA>        NA
1370      <NA>        NA
1371      <NA>        NA
1372      <NA>        NA
1373      <NA>        NA
1374      <NA>        NA
1375      <NA>        NA
1376      <NA>        NA
1377      <NA>        NA
1378      <NA>        NA
1379      <NA>        NA
1380      <NA>        NA
1381      <NA>        NA
1382      <NA>        NA
1383      <NA>        NA
1384      <NA>        NA
1385      <NA>        NA
1386      <NA> 65.698260
1387      <NA> 65.575614
1388      <NA> 65.369538
1389      <NA> 65.480852
1390      <NA> 65.612036
1391      <NA> 65.533773
1392      <NA> 65.513310
1393      <NA> 65.183635
1394      <NA> 65.632586
1395      <NA> 65.431400
1396      <NA> 65.294759
1397      <NA> 65.360027
1398      <NA> 65.328470
1399      <NA> 65.093539
1400      <NA> 65.581364
1401      <NA> 65.367900
1402      <NA> 65.380899
1403      <NA> 65.472944
1404      <NA> 65.388660
1405      <NA> 65.346695
1406      <NA>        NA

plot full predictions

  • What do you think if the predictions? Any better than model 1?
  • Do you agree with the Hmax values generated by the model?
  • If not, why do you think the model is generating those values? (also see next slide)
fit2.predictions %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop), alpha=.2) +
  geom_line(aes(y=prediction, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2", guide="none") +
  facet_wrap(~pop)

DPR data + all predictions

  • This may be helpful in answering the last question
fit2.predictions %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop),
            alpha=.2,
            data = {fit2.predictions %>% filter(pop=="DPR")}) +
  geom_line(aes(y=prediction, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2") 

Let’s try a new model

  • Having Hmax alone vary by population did not produce a very good fit.
  • Exercise: Create and fit a model that allows Hmax and k parameters to vary by population
  • Plot your results: Any better?
  • Perform a model comparison among all 3 models using the loo technique shown above: which fits best?

My formula

f3 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + delta ~ 1,
                  Hmax + k ~ pop,
                  nl=TRUE)

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b     popCC                  Hmax       (vectorized)
               (flat)     b    popDPR                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)
               (flat)     b     popCC                     k       (vectorized)
               (flat)     b    popDPR                     k       (vectorized)

Set the priors

What priors should we use?

prior3 <- c(set_prior("normal(5,3)", lb=0, nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("normal(0, 1)", nlpar = "k", coef="popCC"),
            set_prior("normal(0, 1)", nlpar = "k", coef= "popDPR"),
            set_prior("gamma(20, 3)", nlpar="delta")

)

I use a gamma on k to keep the Intercept k positive.

I use a zero-centered normal for the popCC and popDPR k priors so that these can vary above or below the intercept value of k

model summary

Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         delta ~ 1
         Hmax ~ pop
         k ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      4.04      0.52     3.00     5.04 1.00     3304     2520
delta_Intercept     6.69      0.52     5.73     7.78 1.00     2394     2914
Hmax_Intercept     64.88      1.32    62.42    67.56 1.00     2809     2631
Hmax_popCC        -15.04     21.38   -46.58    30.09 1.00     1897     2174
Hmax_popDPR       -28.29      7.65   -38.87    -8.77 1.00     1738     1277
k_Intercept         0.67      0.01     0.65     0.69 1.00     2614     2394
k_popCC             0.12      0.11    -0.03     0.36 1.00     1702     1690
k_popDPR            0.32      0.07     0.16     0.45 1.00     1644     1414

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.18      0.26    10.69    11.70 1.00     3917     2676

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

model plot

pairs

You can use a pairs plot to see how to parameters co-vary during the fit. Generally we don’t want to see a strong correlation among parameters. Here we see that Hmax and K for popCC have an inverse correlation. This isn’t ideal but I think it is because we have so little CC data. We could impose narrower priors on these parameters but given the lack of data I think I will just leave it.

Predictions

Model comparison

fit1 <- add_criterion(fit1, "loo")
fit2 <- add_criterion(fit2, "loo")
fit3 <- add_criterion(fit3, "loo")
loo_compare(fit1, fit2, fit3)
     elpd_diff se_diff
fit3   0.0       0.0  
fit2  -8.0       3.9  
fit1 -18.1       6.8  

Next

For Monday March 4th

  • Fit a model where delta and Hmax can vary by pop.
  • Fit a model where delta, k, and Hmax can vary by pop. Don’t worry if this doesn’t work well (but if it doesn’t, think about why it doesn’t)
  • Do model diagnostics, plots, and predictive plots on the above.
  • Do a model comparison between all 5 models

Other stuff that Julin is thinking about for future meetings

  • Uncertainty on the predictions
  • On to random effects and all pops
  • CENSORED DATA?